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FOREWORD 


A  literature  search  was  performed  to  assess  the  state  of  the  art  of  the  effort 
to  create  design  criteria  dealing  with  dynamic  buckling  of  cylindrical  shells.  / 
This  work  was  initiated  to  fill  the  need  of  the  Navy  to  have  design  tools  to  use 
in  making  missiles,  mines,  and  sonobuoys.  After  a  comprehensive  search,  articles 
were  critically  reviewed  and  evaluated  on  their  ability  to  produce  design  criteria 
to  meet  the  Navy's  requirements.  Several  recent  papers  indicated  that  the  / 
necessary  theories  were  available  but  had  not  been  applied  to 
parametric  studies  needed. 


producing  the/ 
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INTRODUCTION 


Numerous  ordnance  items,  such  as  missiles,  mines,  and  sonobuoys,  are  subjected 
to  axial  impact  loads  resulting  from  water  entry  or  through-ice  penetration. 

The  use  of  modern  lightweight,  high  strength  materials  for  these  applications 
tends  to  decrease  shell  stability  unless  adequate  stiffening  is  employed.  Accurate 
prediction  of  the  transient  behavior  of  shells  subjected  to  axial  and  side  pulse 
loads  therefore  becomes  of  paramount  importance  if  structural  weight  is  to  be 
minimized. 

At  present  there  are  a  number  of  computer  codes  (STAGS,  for  example)  which 
handle  the  dynamic  buckling  problem  by  discretizing  the  structure  through  the 
use  of  finite  differences  or  elements.  These  codes  can  hardly  be  considered  design 
tools,  because  the  user  must  have  some  idea  of  a  design  prior  to  establishing 
a  model.  In  addition,  these  computer  codes  are  expensive  to  use  and  provide  no 
means  for  cost-effectively  conducting  any  type  of  parametric  design  study;  that 
is,  if  a  given  design  proves  to  be  over-  or  under-designed,  there  is  no  clear 
indication  of  the  options  available  to  the  designer  in  order  to  weight-ef f iciently 
redesign  the  structure. 

A  literature  survey  was  therefore  conducted  to  ascertain  the  current  theo¬ 
retical  methods  available  for  predicting  dynamic  buckling  behavior, 
in  three  parts:  l)a  historical  review,  giving,  in  brief,  a  look  at 
of-the-art,  2)  an  evaluation  of  the  state-of-the-art  with  a  critical 
specific  problem  areas  and  3)  recommendations  for  further  activities 


It  is  presented 
the  state- 
look  at 
in  the  field. 
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HISTORICAL  REVIEW 


The  subject  of  dynamic  buckling  of  cylindrical  shells  has  interested  many 
researchers  for  varying  reasons.  In  the  past  twenty  years  the  aero-space  indus¬ 
try  has  shown  the  greatest  interest  because  of  the  economic  advantages  associated 
with  improved  designs.  The  defense  industry  has  supported  efforts  conmensurate 
with  its  needs  in  missile  and  ordnance  technology.  Others  in  the  scientific 
community  have  received  aid  from  various  organizations  enabling  them  to  advance 
the  state-of-the-art.  Whatever  the  motives,  the  results  include  papers  which 
form  the  foundation  for  further  efforts. 

✓ 

The  work  compiled  herein  can  be  grouped  into  general  categories  in  several 
ways.  One  method  is  by  the  type  of  loading  the  shell  is  subjected  to,  static 
or  dynamic.  Other  ways  include  division  by  type  of  solution  or  by  type  of  shell. 
When  organized  by  type  of  loading  and  put  in  chronological  order,  one  can  trace 
the  general  philosophies  used  in  dealing  with  the  phenomenon.  Any  survey  of  dynam¬ 
ic  buckling  literature  must  contain  literature  of  the  static  case  because  many 
authors  either  simply  extend  static  theory  to  include  the  dynamic  problem,  or 
use  insight  gained  from  the  static  buckling  phenomenon  to  begin  their  study  of 
dynamic  buckling. 

The  following  historical  synopsis  is  organized  by  loading  criteria  and  covers 
literature  published  since  1960.  The  year  1960  has  been  selected  because  of  the 
significant  efforts  begun  in  the  60' s  which  have  become  the  basis  for  present 
day  programs.  The  papers  included  constitute  the  most  significant  programs  of 
the  period.  A  summary  of  literature  sources  is  given  in  Appendix  A.  Notes  taken 
on  each  article  highlighting  Che  important  aspects  of  the  performed  research  are 
presented  in  Appendix  B. 

STATIC  BUCKLING 


Early  researchers  such  as  Donnell,  von  Karman,  Tsien,  and  Koiter  had  developed 
theories,  both  linear  and  non-linear,  concerning  the  buckling  of  circular  cylin¬ 
drical  shells.  These  theories  have  their  basis  in  continuum  mechanics  and  stability. 
Their  strength  lies  in  their  general  form  and  versatility,  allowing  the  inclusion 
of  varying  shell  geometries  and  loading  sequences. 

One  of  the  first  papers  to  be  published  in  the  field  of  static  buckling  of 
shells  in  the  1960's  was  written  by  Brush  and  Morton  in  1961. 1  The  work  was  based 
on  experiment  and  treated  the  shell  as  a  column  with  an  eccentric  load.  The  results 
were  not  important  to  this  study  but  the  simplifying  assumptions  made  give  a  good 
indication  of  the  abilities  of  the  researchers  to  attack  the  problem  in  a  realistic 
sumner  at  that  time. 


*Brush,  D.  0.  and  Morton,  F.  G.,  "Stress  &  Stability  Analysis  of  Cylindrical  and 
Conical  Shells,"  Tech.  Report  LMSD-894808,  Lockheed  Missiles  and  Space  Div., 
Sunnyvale,  CA,  31  May  1961. 
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In  attempting  to  analyze  the  complex  combined  forces  of  axial  compression 
and  external  pressure  encountered  in  launch  vehicle  structures,  Lakshmikantham 
and  Gerard,  in  1964,  wrote  an  analysis  of  unstiffened  circular  cylindrical  shells 
based  upon  the  Donnell  equation  in  terms  of  normal  displacements . 2  The  analysis 
included  initial  imperfections  and  viewed  the  buckling  stress  in  terms  of  the 
Batdorf  parameter,  Z.  Imperfections  or  prebuckling  deformations  were  noted  to 
have  an  unusual  influence  on  the  structures,  particularly  when  an  axial  load  was 
combined  with  at  least  a  small  circumferential  load,  i.e.,  axial  loads  made  the 
shell  imperfection  sensitive  while  the  shell  was  relatively  insensitive  to  imper¬ 
fections  under  a  pure  circumferential  load. 

The  method  of  Singer,  Baruch,  and  Harari  which  would  be  used  by  others  later, 
dealt  with  shell  ring  and  stringer  stiffeners. 3  It  was  the  first  to  "smear"  the 
equivalent  stiffening  effect  of  the  stiffener  over  the  entire  shell.  It  is  inter¬ 
esting  that  for  purposes  of  analysis  it  was  assumed  that  stiffeners  were  closely 
spaced  so  that  the  smearing  took  care  of  both  inside  and  outside  effects,  even 
though  the  author  concluded  that  there  were  considerable  differences  between  inside 
and  outside  stiffeners.  Small  deflection  theory  using  Donnell  type  stability 
equations  was  utilized. 

In  1967  the  direct  method  of  Liapunov  using  Donnell-type  shell  theory  was 
used  to  analyze  a  specific  shell  by  making  many  simplifying  assumptions.2*  But 
it  was  Hutchinson  and  Amazigo^  who  made  advances  using  non-linear  Donnell-type 
strain  displacement  equations  analyzing  imperfection  sensitivity  with  respect 
to  Z  for  stiffened  shells.  The  solution  could  not  be  applied  to  design  a  shell 
but,  again,  much  insight  was  gained.  One  result  of  their  investigation  was  the 
prediction  of  a  particular  imperfection  sensitivity  in  the  low  Z  range''/ 

The  next  significant  paper  was  written  by  Horton  and  Tenerelli.6  They 
presented  both  analytical  and  experimental  work  dealing  with  ^imperfection  sensi¬ 
tivity  under  axial  and  compressive  radial  loads.  Experimentally  they  showed  that 

2  ■  .  y  . 

Lakshmikantham,  C.  and  Gerard,  G.,  "Elastic  Stability  of  Cylindrical  Shells  Under 
Axial  and  Lateral  Loads,"  Journal  of  Royal  Aerosp^e  Society  Vol.  68,  Nov  1964, 
pp.  773-775.  / 

/ 

/ 

3 

Singer,  J.,  Baruch,  M.  and  Harari,  0.,  "Op^the  Stability  of  Eccentrically 
Stiffened  Cylindrical  Shells  Under  Axiqf  Compression,"  International  Journal  of 
Solids  and  Structures,  Vol.  3,  Oct  196^,  pp.  445-470. 

^Hegemier,  G.  A.,  "Stability  of  Cylindrical  Shells  Under  Moving  Loads  by  the 
Direct  Method  of  Liapunov,"  Office  of  Naval  Research,  Contract  N00014-67-A-0109- 
0003,  Task  NR  064-496  Tech./fteport  No.  5,  Dept,  of  Aerospace  and  Mech.  Eng.  Sci., 
Univ.  Cal.,  San  Diego,  Jvt  1967. 

^Hutchinson,  J.  W.  aptd  Amazigo,  J.  C.,  "Imperfection  Sensitivity  of  Eccentrically 
Stiffened  Cylindrical  Shells,"  AIAA  Journal,  Vol.  5,  March  1967,  pp.  392-401. 

/ 

^Horton,  W.  H-and  Tenerelii,  0.  J.,  "The  Instability  of  Unstiffened  and  Ring 
Stabilized /thin  Walled  Circular  Cylinders  Under  Non-Uniform  Axial  Load 
Conditions,"  Israel  Journal  of  Technology,  Vol.  6,  No.  1-2,  1968,  pp.  11-24. 
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Che  critical  compressive  stress  was  30  to  402  less  than  the  critical  bending 
stress  and  that  due  to  imperfections,  an  increased  allowable  stress  could  be 
used  in  flexure  if  it  were  possible  to  assess  the  probability  that  the  location 
of  the  weakest  section  of  the  shell  would  not  correspond  to  the  region  of  maximum 
stress.  The  analysis  was  accomplished  using  a  computer  program  which  solved 
the  three  simultaneous  equations  of  Flugge. 

H.  A.  Shumik  published  a  paper^  claiming  to  attack  the  dynamic  stability  of 
cylindrical  shells  under  the  action  of  dynamic  radial  pressure.  He  assumed  inertia 
could  be  neglected  thereby  making  it  a  quasi-static  problem.  This  depicts  a 
common  problem  associated  with  buckling  analysis,  that  is,  the  results  of 
assumptions.  El  Naschie^  encountered  this  same  hurdle  in  his  mathematical  treatment 
of  buckling  patterns  using  a  two  dimensional  Riemann  space  model. 

The  current  state  of  the  art  in  static  buckling  analysis  is  presented  by  Weller 
and  Singer. ^  By  comparing  experimental  and  theoretical  results  it  is  shown  that 
linear  theory  can  be  used  effectively  as  a  first  approximation  for  predicting 
buckling  loads  of  stringer  stiffened  shells.  For  the  analysis,  general  instability 
was  calculated  using  the  concept  of  "smeared"  stiffeners.  The  results  of  the 
experimentation  refuted  the  prediction  by  Hutchinson  of  an  imperfection 
sensitivity  in  the  low  Z  range. 

DYNAMIC  BUCKLING 

One  of  the  first  American  papers  to  appear  in  the  field  of  buckling  of 
cylindrical  shells  under  dynamic  loads  was  written  by  J.  C.  Yao  in  1961. 10  This 
work  was  followed  by  another  in  1963. ^  Both  were  analytical  using  Mathieu's 
equation  and  a  modified  Donnell  equation  of  motion  to  examine  the  stability  of  the 
system  under  certain  assumptions  by  seeking  two  wave  parameters.  For  his  first 
paper  the  flexural  mode  was  coupled  with  the  axisymmetric  mode  and  was  shown  to 
respond  in  an  oscillatory  or  exponential  fashion  depending  upon  the  magnitude 
of  the  square  wave  pressure  pulse.  In  the  second  analysis  periodic  load 
frequencies  were  assumed  to  be  very  much  less  than  the  frequencies  of  the  stress 

^Shumik,  M.  A.,  "The  Stability  of  Cylindrical  Shells  Under  the  Action  of  Dynamic 
Radial  Pressure,"  Foreign  Technology  Division,  Wright-Patterson  AFB,  Ohio, 

25  Jun  1974,  13  pp. 

8 

El  Naschie,  M.  S.,  "Localized  Diamond-Shaped  Buckling  Patterns  of  Axially 
Compressed  Cylindrical  Shells,"  AIAA  Journal,  13,  6,  pp.  837-838,  (Technical  Notes), 
Jun  1975. 

o 

Weller,  T.  and  Singer,  J.,  "Experimental  Studies  on  the  Buckling  Under  Axial 
Compression  of  Integrally  Stringer  Stiffened  Circular  Cylindrical  Shells," 

Journal  of  Applied  Mechanics,  44,  Dec  1977,  pp.  721-730. 

*®Yao,  John  C.,  "Stability  of  a  Cylinder  Under  Dynamic  Radial  Pressure,"  A.R.S. 
Journal,  Vol.  31,  No.  12,  Dec  1961,  pp.  1705-1708. 

**Yao,  John  C. ,  "Dynamic  Stability  of  Cylindrical  Shells  Under  Static  and  Periodic 
Axial  and  Radial  Loads,"  AIAA  Journal,  Vol.  1,  No.  6,  Jun  1963,  pp.  1391-1396. 


NSWC  TR  79-447 


waves  in  the  cylinder  and  the  cylinders  were  of  medium  length,  (10^  <  Z  <  10^). 
Regions  of  stability  were  examined  for  cylinders  with  no  stiffeners  where  the 
form  of  the  displacement  was  assumed.  End  conditions  were  ignored  at  this  point 
whereas  by  1977  they  would  be  considered  important  parameters. 9 

The  question  of  whether  to  consider  the  effects  of  inertia  in  the  dynamic 
problem  was  generally  answered  by  making  assumptions.  Roth  and  Klosner  in  1964 1 
assumed  the  effects  of  inertia  negligible  in  their  study  of  the  nonlinear  response 
of  cylindrical  shells  subjected  to  axial  loads  of  constant  magnitudes  of  prolonged 
or  finite  time  duration.  The  analysis  used  Donnell's  strain-displacement  relation¬ 
ships,  Hooke's  law  for  thin  shells,  equations  for  elastic,  extensional  strain 
energy,  Hamilton's  principle,  and  an  Airy  stress  function  to  give  two  nonlinear 
differential  equations  in  terms  of  the  stress  function  and  radial  displacement. 

This  method  resulted  in  the  first  data  related  to  impulse  loading.  The  critical 
impulse  was  found  to  decrease  as  the  time  duration  of  loading  decreased  and  the 
dynamic  buckling  stress  decreased  with  an  increase  in  the  time  duration  of  the 
loading.  The  lack  of  inertia  terms,  however,  leave  some  room  for  doubt  concerning 
the  accuracy  of  the  analysis  for  all  loads  and  shells. 

The  first  significant  paper  concerning  cylindrical  shells  under  axial  impact 
which  includes  experimentation  was  written  by  Linberg  and  Herbert  in  1966.^ 

Linear  theory  was  applied  to  allow  the  buckling  form  to  assume  its  own  shape. 
Simplified  Donnell-type  equations  were  used  which  incorporated  imperfections. 

The  experimentation  was  run  simultaneously  with  the  theoretical  investigation 
to  see  if  the  shell,  once  deformed  into  its  initial  deflection  pattern  defined 
by  linear  shell  theory,  would  continue  in  that  form.  The  experimentation  used 
a  unique  method  of  applying  explosives  to  accelerate  the  shell  while  a  high  speed 
camera  monitored  its  motion  at  240000  frames  per  second.  The  loading  was  uniform 
and  controlled  to  within  2  microseconds.  The  experimentation  confirmed  the  pre¬ 
mise  that  early  small-deflection  buckling  dictates  the  pattern  into  which  large 
deflection  buckles  form. 

14 

Further  experimentation  was  performed  by  Coppa  using  cylinders  set  on  a 
stand  pressurized  with  a  free  end  and  impacted  by  a  drop  test  apparatus.  The 
test  cylinders  were  unstiffened  and  assumed  perfect.  Like  the  article  by  Roth 
and  Klosnerl2  the  initial  mode  of  buckling  was  examined  and  certain  insights  were 


See  footnote  9  on  page  10. 

12 

Roth,  R.  S.  and  Klosner,  J.  M. ,  "Nonlinear  Response  of  Cylindrical  Shells 
Subjected  to  Dynamic  Axial  Load,"  AIAA  Journal,  Vol.  2,  Oct  1964,  pp.  1788-1794. 

13 

Linberg,  H.  and  Herbert,  R.  E.,  "Dynamic  Buckling  of  Thin  Cylindrical  Shells 
Under  Axial  Impact,"  Journal  of  Applied  Mechanics,  Vol.  33,  Jan  1966,  pp.  105-112. 

14 

Coppa,  A.  P. ,  "Effect  of  End  Conditions  on  the  Buckling  of  Cylindrical  Shells 
Under  Axial  Compressive  Impact,"  G.E.  Missile  and  Space  Div.,  Space  Sciences 
Laboratory,  R66SD24,  May  1966,  22  pp. 


11 


^  :  «*»>•  L-*r. 
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gained.  It  was  postulated  that  triangular  buckles  represent  a  lower  energy  mode 
of  response  than  the  symmetric  mode.  Symmetrical  buckles  result  from  predominantly 
outward  displacements,  triangular  buckles  inward.  When  a  shell  is  impacted,  it 
is  forced  to  respond  quickly  in  an  energetically  efficient  deformation  mode. 
Therefore  the  dominant  mode  is  symmetric  initially.  Whether  this  mode  prevails 
depends  upon:  shell  thickness,  amount  of  lateral  support,  and  the  relationship 
between  impacted  strain  and  the  yield  strain  of  the  material.  It  was  also  noted 
that  triangular  buckling  was  independent  of  end  conditions  which  had  only  been 
hypothesized  previously. 

The  first  paper  to  include  inertia,  imperfections,  and  stiffeners,  dealing 
with  dynamic  buckling  used  large  deflection  theory  and  assumed  clamped  ends.l* 

The  "smearing"  technique  was  used  to  deal  with  the  stiffeners  and  assumed  forms 
were  made  for  radial  displacements  and  end  shortening. 

Under  the  guidance  of  Professer  R.  E.  Ball,  C.  C.  Cromer^  considered  the 
nonlinear  dynamic  response  of  ring  stiffened  shells.  Inertia  was  not  included 
but  imperfections  were  in  the  nonlinear  Donnell-type  equations.  Damping  due  to 
viscous  forces  was  also  incorporated.  This  resulted  in  a  lengthy  computer  program 
for  ring  stiffened  shells  under  axisymmetr ic  radial  loads  only. 

Another  computer  program  using  a  modified  membrane  theory  including  transverse 
shear  deformation  was  written  by  Mortimer,  Rose,  and  Chou  in  1969.  ^  Tests  were 
also  conducted  to  verify  the  computer  output  using  a  pendulum  impactor,  photo¬ 
diodes  and  light  source.  The  analysis  includes  no  inertia  or  imperfection  effects 
and  was  for  simple  unstiffened  shells. 

1 8 

In  1972,  Professer  R.  E.  Ball  presented  a  computer  program  using  Sander's  non¬ 
linear  thin  shell  theory  for  conditions  of  small  strains  and  moderately  small  rota¬ 
tions.  The  program  could  analyze  shells  of  revolution  for  any  arbitrary  load 
definable  by  a  Fourier  series.  It  approximated  dependent  variables  using  Fourier 
series  and  treated  nonlinear  modal  coupling  with  pseudo  loads.  Limitations  included 
its  application  only  when  geometric  and  material  properties  were  axisymmetric. 


Dietz,  W.  K.,  "On  the  Dynamic  Stability  of  Eccentrically  Reinforced  Circular 
Cylindrical  Shells  Technical  Report,  Sep  1,  1965  -  Jan  31,  1967,"  Syracuse 
University  Research  Institute,  NY,  Applied  Mechanics  Lab.,  31  Jan  67,  190  pp. 


Cromer,  Charles  Clinton,  "An  Investigation  of  the  Nonlinear  Dynamic  Response  of 
Cylindrical  Shells  Under  Transient  Pressure,"  Naval  Post  Graduate  School 
Monterey,  CA,  April  1969,  281  pp.  ’ 

t 

Chou,  P.  C. ,  Mortimer,  R.  W.  and  Rose,  J.  L. ,  "Wave  Propagation  in  Stepped  and 
Joined  Shells  Annual  Report,  1  Sep  1969,"  Drexel  Institute  of  Technology, 
Philadelpha,  PA,  Wave  Propagation  Research  Center,  24  Oct  1969,  37  pp. 

Ball,  R.  E.,  "A  Program  for  the  Nonlinear  Static  and  Dynamic  Analysis  of 
Arbitrarily  Loaded  Shells  of  Revolution,"  Computers  and  Structures,  Vol.  2, 
pp.  141-162,  Pergamon  Press  1972. 
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the  material  was  isotropic,  and  shells  had  no  stringers.  The  results  demanded 
interpretation  to  ascertain  whether  the  shell  reached  its  buckling  mode.  Subse¬ 
quent  articles  by  Professer  Ball  dealt  with  similar  subject  matter  but  were  not 
pertinent  to  this  study . 19, 20, 21 

Other  analytic  approaches  included  the  use  of  extremely  general  field 
equations^  and  an  energy  approach  to  examine  the  stability  of  arches  and  spherical 
caps. 23  jn  the  analysis  using  the  energy  method  it  was  mentioned  that  complex 
structures,  imperfections,  and  real  material  properties  pose  problems  to  the 
successful  use  of  this  approach. 

The  4  degree  of  freedom  model,  the  nonlinear  Donnell  equations,  and  the  basic 
method  of  solution  performed  by  Roth  and  Klosner^  in  1964  were  reworked  by 
Y.  S.  Tamura  in  1973.24  As  with  Roth  and  Klosner,  stiffeners  were  neglected. 
However,  unlike  them,  inertia  was  approximated,  damping  of  the  axial  motion  was 
included,  and  parametric  studies  of  wave  number,  mass  of  loaded  edge,  and 
damping  factor  were  conducted.  It  was  concluded  that  due  to  frequency  coupling 
between  axial  and  radial  motions,  axial  inertia  plays  an  essential  role  in  charac¬ 
terizing  the  dynamic  instability  of  a  finite  length  shell.  This  was  in  direct 
opposition  to  all  previous  assumptions  about  inertia  effects. 


Ball,  R.  E.,  "Survivability  of  the  Five-Inch  Gun  Launched  Finned  Motor  Case," 
Naval  Post  Graduate  School,  Monterey,  CA  NPS-57Bp72081A,  10  August  1972. 

20 

Ball,  R.  E.  and  Burt,  J.  A.,  "Dynamic  Buckling  of  Shallow  Spherical  Shells," 

J.  Applied  Mechanics,  June  1973,  pp.  411-416. 

21 

Ball,  R.  E. ,  et  al.,  "A  Comparison  of  Computer  Results  for  the  Dynamic  Response 
of  the  LMSC  Truncated  Cone,"  Computers  and  Structures,  Vol.  4,  pp.  485-498, 
Pergamon  Press  1974. 

22 

Darevskii,  V.  M.  ,  "Stability  of  a  Cylindrical  Shell  to  Dynamic  Axial  Loading," 
Mekh.  Tverd,  Tela,  (USSR),  No.  2,  1973,  pp.  162-172. 

23 

Bhatia,  P.  and  Babcock,  C.  D. ,  Jr.,  "Dynamic  Buckling  of  Structures," 

California  Institute  of  Technology,  Pasadena  Graduate  Aerospace  Labs.  Feb  1973, 
64  pp. 

12 

See  footnote  12  on  page  11. 

24 

Tamura,  Y.  S.,  "Dynamic  Stability  of  Cylindrical  Shells  Under  Step  Loading," 
Air  Force  Office  of  Scientific  Research,  Calif.  Inst,  of  Tech.,  AD-777-291, 

June  1973. 


13 


NSWC  TR  79-447 


Approaching  the  subject  from  an  elastic  systems  viewpoint  the  nonlinear  von 
Karman-Donnell  equations  can  lead  to  a  general  analysis  of  the  phenomenon  of 
resonance  and  its  importance  to  system  stability.  One  such  analysis  was  conducted 
by  Lange^S  of  a  specific  shell  structure  to  show  that  nonlinear  coupling  of  normal 
modes  of  the  solutions  to  the  equilibrium  equations  can  have  significant  bearing 
on  the  stability  of  elastic  systems. 

26 

C.  A.  Fisher  and  C.  W.  Bert  developed  a  computer  program  which  neglected 
imperfections  and  treated  outside  stiffeners  perfunctorily,  but  included  discrete 
rings  and  stiffeners  subjected  to  rapidly  applied  axial  compressive  loads.  Both 
general  and  localized  instability  were  considered.  The  major  drawback  was  that 
the  program  used  a  tremendous  amount  of  computer  time  and  core  space. 

24 

Publishing  at  the  same  time  as  Tamura  ,  Tsui  and  Lakshmikantham  had  written 
their  first  in  a  series  of  articles  dealing  with  the  nonlinear  and  linear  response 
of  stiffened  shells  to  dynamic  loads. 27  Tamura  had  worked  for  the  Air  Force 
Office  of  Scientific  Research  and  Tsui  and  Lakshmikantham  for  the  Army  Materials 
and  Mechanics  Research  Center.  Independently  they  used  the  same  approach  to  solve 
very  similar  problems:  Tamura,  cylindrical  shells;  Tsui,  stiffened  shells.  The 
articles  by  Tsui  and  Lakshmikantham  investigated  the  use  of  the  von  Karman-Donnell 
non-linear  approach27]  the  linear  theory  of  Singer^®,  and  again  the  Donnell  type 
formulation . 29 , 30  their  first  paper  they  indicated  that  dynamic  buckling  loads 

are  less  for  inside  stiffened  shells,  a  result  similar  to  that  of  Hutchinson  and 

25Lange,  C.  G.  and  Newell,  A.  C.,  "The  Role  of  Resonant  Interactions  in  the  Dynamic 
Stability  of  Elastic  Systems,"  California  University,  Los  Angeles,  Department  of 
Mathematics,  1973,  13  pp. 


9 

Fisher,  C.  A.  and  Bert,  C.  W. ,  "Design  of  Crashworthy  Aircraft  Cabins  Based  on 
Dynamic  Buckling,"  J.  Aircraft,  Vol.  10,  No.  11,  November  1973,  pp.  693-695. 


See  footnote  24  on  page  13. 

27  ... 

Tsui  and  Lakshmikantham,  "The  Nonlinear  Dynamic  Response  of  Stiffened  Shells 

Under  Compressive  Axial  Impulse,"  Applied  Mechanics,  4th  Canadian  Congress, 
(CANCAM  73)  2  May  -  1  June  1973,  Montreal,  Canada. 

28 

Tieu-Yu  Tsui,  et  al,  "Linear  Stability  of  Eccentrically  Stiffened  Cylindrical 
Shells  Under  Axial  Compression,"  Army  Mat.  &  Mech.  Research  Center,  Watertown, 
Mass.,  AD-775-445,  Dec.  1973. 


29  , 

Lakshmikantham,  C.  and  Tsui,  T.,  "Dynamic  Stability  of  Axially-Stiffened 

Imperfect  Cylindrical  Shells  Under  Axial  Step  Loading,"  AIAA  Journal,  Vol.  12, 

Feb  1974,  pp.  163-169. 


30  ..... 

Lakshmikanthan,  C.  and  Tsui,  T. ,  "Dynamic  Buckling  of  Ring  Stiffened  Cylindrical 

Shells,"  Army  Materials  and  Mech.  Research  Center,  Watertown,  MA,  AIAA 

Journal,  Vol.  13,  No.  9,  Sep  1975,  pp.  1165-70. 
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Amaz igc  for  the  static  case.  It  is  in  these  papers  that  is  found  the  first 
efforts  to  quantify  the  effects  of  stiffener  eccentricity.  Eccentricity  is  defined 
as  the  measure  of  the  distance  between  the  centroid  of  the  stiffener  and  the 
middle  surface  of  the  shell  with  sign  (+  or  -)  indicating  inside  or  outside  orienta¬ 
tion.  In  all  cases  the  computer  was  used  to  do  parametric  studies  using  tremendous 
amounts  of  computer  time. 

30 

Just  prior  to  the  last  of  the  articles  by  Tsui  and  Lakshmikantham  Tamura 
co-authored  a  paper  with  Babcock-^  dealing  again  with  the  subject  of  dynamic  buckling 
of  unstiffened  cylindrical  shells.  Both  treatises30,31  dealt  with  the  subject 
in  almost  identical  fashions:  using  nonlinear  Donnell's  equations  in  terms  of 
displacements,  assuming  radial  displacement  and  initial  imperfection  functions, 
they  evaluated  the  energy  of  the  system,  used  Hamilton's  principle,  and  then  used 
an  approximation  method  to  integrate  the  resulting  nonlinear  ordinary  differential 
equations.  In  both  cases  axial  inertia  was  found  to  have  an  important  influence 
on  the  final  buckling  characteristics  of  the  shell.  Tamura  also  performed  a  static 
analysis  and,  when  comparing  it  to  the  dynamic  analysis,  found  that  for  finite 
length  shells  reflections  of  the  axial  stress  wave  must  be  taken  into  account. 

New  innovations  in  the  Tamura  paper  included  assessment  of  the  influence  of 
a  realistic  triggering  mechanism  on  the  calculated  dynamic  buckling  loads  and 
the  inclusion  of  in-plane  inertia.  Together  with  the  progress  in  the  areas  of 
axial  inertia,  imperfection  approximations,  application  of  Donnell's  equations, 
and  numerical  integration  techniques  this  represents  the  most  thorough  nonlinear 
analysis  of  unstiffened  cylindrical  shells  under  a  dynamic  load  to  date. 

Lakshmikantham  and  Tsui  have  performed  the  most  comprehensive  nonlinear  analysis 
of  stiffened  cylindrical  shells  under  a  dynamic  load  to  date. 30 

The  question  of  critical  impact  velocity  for  a  cylindrical  shell  was  raised 
by  A.  B.  Efimov. 32  He  derived  a  formula  which  calculated  the  critical  velocity 
of  a  cylindrical  shell  on  longitudinal  impact  against  an  infinite  mass.  The  theory 
was  based  upon  a  nonlinear  von  Karman-Donnell  equation  in  terms  of  radial  displace¬ 
ment.  The  author  claimed  that  the  formula  agreed  with  experimentation  to  within 
12%.  No  proof  of  experimentation  was  presented.  Efimov's  mentor,  Grigolyuk, 


^See  footnote  5  on  page  9. 

30 

Lakshmikanthan,  C.  and  Tsui,  T.,  "Dynamic  Buckling  of  Ring  Stiffened  Cylindrical 
Shells,"  Army  Materials  and  Mech.  Research  Center,  Watertown,  MA,  AIAA  Journal, 

Vol .  13,  No.  9,  Sep  1975,  pp.  1165-70. 

31 

Tamura,  Y.  S.  and  Babcock,  C.  D.,  "Dynamic  Stability  of  Cylindrical  Shells  Under 
Step  Loading,"  Transactions  of  the  ASME,  March  1975,  pp.  190-194. 

32 

Efimov,  A.  B.,  Malyi,  V.  I.  and  Uteshev,  S.  A.,  "Loss  of  Stability  of  a  Cylindrical 
Shell  on  Longitudinal  Impact,"  Foreign  Technology  Division,  Wright-Patterson 
AFB,  Ohio,  27  Jun  1974,  12  pp. 
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published  a  paper  in  1975  which  was  written  in  the  same  style,  not  displaying 
any  evidence  of  real  work,  but  investigating  the  effect  of  external  pressure 
impulse  shape  on  circular  cylindrical  shells. 33 

Three  studies  followed  using  linear  buckling  theory. The  first  in 
1974  considered  arbitrary  axisymmetric  shells  subject  to  axisymmetric  periodic 
loads  of  long  duration,  thereby  creating  a  pseudo-dynamic  loading  sequence. 34 
The  second  in  1975  was  filled  with  assumptions  to  disregard  higher  order  terms 
to  lead  to  Mathieu's  equation  in  order  to  generate  parametric  curves  showing 
minimum  buckling  loads. 35  The  last,  by  Maymon  and  Libai36  in  1977,  linearized 
the  Donnell-type  equations  derived  by  Baruch  and  Singer  to  yield  a  computer  program 
to  analyze  a  shell  under  axial  impact.  Similar  to  others  who  have  dealt  with 
stiffened  shells,  a  smearing  technique  was  used.  The  results  included  a  parametric 
study  of  different  shells  and  an  analysis  of  directions  to  take  in  further  work. 

36 

The  paper  of  Maymon  and  Libai  was  quickly  followd  by  that  of  Libai  and 
Durban. 3'  Instead  of  following  his  own  recommendations,  Libai 's  second  paper 
concerned  itself  with  circumferentially  varying  axial  edge  loads  or  thermal  loads 
as  opposed  to  axial  impact.  The  problem  was  also  approached  differently  than 
before.  It  was  treated  as  an  eigenvalue  problem  derived  from  the  linearized 
Donnell  equations  and  membrane  theory.  If  Libai  had  only  continued  his  efforts 
along  the  lines  of  axial  impact  perhaps  design  criteria  for  shells  subject  to 
this  loading  could  have  been  established. 

Perhaps  the  most  comprehensive  treatise  of  all  those  found  in  the  literature 
thus  far  is  that  of  D.  G.  Zimcik.38  His  paper  reports  the  results  of  experimental 
testing  of  thin  walled  circular  cylindrical  shells  subjected  to  dynamic  transient. 


33 

Grigolyuk,  E.  1.  and  Srebovski,  A.  I.,  "Effect  of  Shape  of  External  Pressure 
Impulse  and  Design  Parameters  of  Thin  Circular  Cylindrical  Shells  on  Their 
Stability,"  NISC,  Washington,  DC,  Translation  Division,  Dec  1975,  8  pp. 

3^ 

xalnins,  A.,  "Dynamic  Buckling  of  Axisymmetric  Shells,"  Trans.  ASME  Ser.  E., 

Vol.  4,  Dec  1974,  pp.  1063-68. 

35 

Goroshko,  D.  A.  and  Emel-Yanenko,  V.  V. ,  "Dynamic  Stability  of  Layered  Anisotropic 
Shells,"  Prikl ,  Mekh.  (USSR),  Vol.  11,  No.  7,  July  1975,  pp.  42-49. 
or:  Sov.  Appl.  Mech.  (USA),  Vol.  11,  No.  7,  720-25,  Jul  1975. 

36 

Maymon,  G.  and  Libai,  A.,  "Dynamics  and  Failure  of  Cylindrical  Shells  Subjected 
to  Axial  Impact,"  AIAA  Journal,  Vol.  15,  No.  11,  Nov  1977,  pp.  1624-1630. 

37 

Libai,  A.  and  Durban,  D. ,  "Buckling  of  Cylindrical  Shells  Subject  to  Nonuniform 
Axial  Loads,"  Journal  of  Applied  Mechanics,  44,  Dec  1977,  pp.  714-720. 

38 

Zimcik,  D.  G. ,  "Stability  of  Circular  Cylindrical  Shells  Under  Transient  Axial 
Impulse  Loading,"  Institute  for  Aerospace  Studies,  Toronto  University, 

Mar  1976,  137  pp. 
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axial,  square  wave  loading  of  varying  time  duration.  Special  manufacturing 
techniques  were  used  to  make  three  types  of  shells  of  birefringent  liquid  epoxy. 
The  types  included  were:  geometrically  near  perfect,  axisymmetric  imperfect, 
and  asymmetric  imperfect.  Shell  response  was  used  to  verify  predictions  obtained 
from  both  linear  and  nonlinear  analysis  based  upon  von  Karman-Donnell  equations 
and  a  modified  linear  stress  wave  analysis.  Experimentation  and  analysis  are 
claimed  to  show  good  agreement. 

In  view  of  the  fact  that  the  objective  of  this  study  is  to  develop  analytical 
methods  to  be  used  for  generating  a  useful  design  tool,  the  literature  search 
was  expanded  to  include  the  Computer  Software  Management  and  Information  Center 
(COSMIC). 

COSMIC  is  operated  for  NASA  by  the  University  of  Georgia  for  the  purpose  of 
making  computer  programs  developed  in  the  space  program  available  to  the  public. 
The  prices  charged  for  programs  by  COSMIC  are  established  in  accordance  with  NASA 
policy  to  recover  as  large  a  portion  of  COSMIC's  operating  expenses  as  possible, 
without  making  programs  prohibitively  expensive  for  small  firms.  In  actual 
practice,  NASA  subsidizes  about  one-third  of  the  cost  of  the  services  provided 
by  COSMIC. 

An  exhaustive  survey  of  the  COSMIC  program  library  found  several  programs  of 
interest.  Abstracts  of  two  dynamic  programs  and  one  static  program  which  hold 
much  promise  for  this  study  are  included  in  Appendix  C. 

LITERATURE  EVALUATION 

While  the  preceeding  literature  survey  showed  that  a  great  deal  of  work  has 
been  done  on  the  subjects  of  static  and  dynamic  buckling,  it  can  be  seen  that 
most  of  the  work  has  been  left  open  ended  or  completed  only  for  the  simpler  of 
shell  wall  geometries.  It  is  no  surprise  that  this  is  a  prevalent  finding  since 
the  subject  is  so  encompassing  that  it  requires  the  resolution  of  complex 
multivariable  relationships  that  usually  takes  years  to  complete. 

Initially  the  major  stumbling  block  for  researchers  was  the  creation  of  a 
mathematical  model  of  sufficient  complexity  and  degrees  of  freedom  to  include 
all  of  the  governing  parameters.  The  work  of  Roth  and  Klosner^  hurdled  this 
obstacle  using  a  Donnell-type  equation,  Hamilton's  principle,  the  introduction 
of  an  Airy  stress  function,  and  subsequent  use  of  a  Runge-Kutta  numeric  technique. 
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A  copy  of  their  paper  is  included  in  Appendix  E.  Most  studies  published  after 
this  one  used  this  technique  while  incorporating  or  modifying  differing 
assumptions .13,16,24,25,27,29,30,31,32,36,38 

Inertia  terms  have  been  a  source  of  conflict.10,11’25*29, 30,31,34,38 
Lakshmikantham  and  Tsui  finally  admitted  that  "normal"  inertia  was  important  enough 
to  include  in  their  final  paper,  long  after  Tamura  and  Babcock  had  shown  its 
importance. 24, 31  Zimcik- ^  included  axial,  rotary  and  radial  inertia  in  his  work. 

It  is  the  opinion  of  these  experts  that  axial  inertia  mist  be  included  in  any 
dynamic  buckling  model. 


footnote 

13 

on 

page 

11. 

16See 

footnote 

16 

on 

page 

12. 

24See 

footnote 

24 

on 

page 

13. 

25See 

footnote 

25 

on 

page 

14. 

27See 

footnote 

27 

on 

page 

14. 

29See 

footnote 

29 

on 

page 

14. 

30See 

footnote 

30 

on 

page 

14. 

31See 

footnote 

31 

on 

page 

15. 

32See 

footnote 

32 

on 

page 

15. 

36See 

footnote 

36 

on 

page 

16. 

38See 

footnote 

38 

on 

page 

16. 

10- 

See 

footnote 

10 

on 

page 

10. 

11  See 

footnote 

11 

on 

page 

10. 

34  See 

footnote 

34 

on 

page 

16. 

18 


NSWC  TR  79-447 


Initial  imperfections  effect  the  buckling  load  and 
mode. 2, 5, 6, 12, 13, 16, 23, 24,27, 29, 30, 31, 32, 36, 38  The  debate  concerning  imperfections 

centers  on  the  method  of  including  them  in  the  analysis,  not  on  whether  they  should 
or  should  not  be  included.  The  best  solutions  were  presented  by  Tamura  and 
Babcock^l  who  used  experimental  data,  and  Zimcik^S  who  made  and  measured  his  own. 

In  general  the  problem  was  dealt  with  by  assuming  some  form  for  the  imperfection. 

A  summary  of  these  methods  is  presented  in  Appendix  D.  One  of  the  results  of 
the  Maymon-Libai  paper^6  presented  the  current  opinion  best  saying,  "more  effort 
is  needed  in  establishing  real  statistical  imperfection  distributions."  Any 
further  effort  in  the  field  must  recognize  this  problem  and  deal  with  it. 
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The  sole  use  of  an  energy  method  has  been  ruled  out  except  in  the  evaluation 
of  instability  modes. ^3  All  current  work  is  being  done  using  the  equations  of 
motion  based  on  Donnell's  equations,  using  energy  methods  only  to  satisfy  boundary 
equations . 

The  results  of  the  best  papers  were  very  similar  to  the  objective  of  this 
study,  that  of  the  generation  of  parametric  curves.  The  similarities,  however, 
exist  only  in  general.  In  all  cases,  the  programs  used  to  perform  the  studies 
by  these  authors  were  either  very  time  consuming,  requiring  large  amounts  of  core 
storage,  or  written  in  a  too  specific  manner,  not  easily  changed  for  the  purposes 
of  dealing  with  the  objective.  In  no  case  had  the  parametric  studies  been  intended 
for  anything  more  than  qualitative  analysis. 

The  exhaustive  survey  of  computer  programs  showed  that  even  COSMIC  does  not 
have  the  type  of  program  needed  to  perform  the  parametric  study.  There  were  a 
few  programs  which  appear  to  lend  themselves  to  being  used  as  "core”  programs 
containing  the  appropriate  theories.  These  could  be  set  up  to  run  with  a  calling 
preprocessor  to  generate  parametric  curves.  However,  as  yet  this  has  not  been 
done . 


RECOMMENDATIONS  FOR  FURTHER  EFFORTS 

Based  upon  the  research  of  previous  authors  presented  here,  further  efforts 
in  this  area  may  be  extremely  fruitful.  The  analytical  ground  work,  in  terms 
of  the  general  equations  of  motion  and  methods  of  solution,  has  been  conducted, 
and  several  computer  programs  are  available  utilizing  this  knowledge.  A  further 
development  effort  could  conceivably  be  to  convert  one  of  these  programs  into 
a  design  program  which  would  vary  a  few  of  the  geometric  parameters,  regarding 
weight  as  a  merit  function  in  the  design  variable  space,  and  then  to  determine 
the  minimum  weight  design  by  mathematical  programming  techniques. 

To  determine  how  these  programs  may  be  made  useful  in  developing  design  and 
trade-off  curves,  a  few  of  the  programs  and  their  associated  documentation  were 
purchased.  Using  the  literature  survey  as  background  on  the  important  aspects 
and  acceptable  theories  of  dynamic  buckling,  it  was  concluded  that  modification 
of  one  or  more  of  these  programs  to  yield  design  codes  is  possible.  The  abstracts 
of  two  dynamic  programs  and  one  static  program  that  hold  good  promise  for  this 
work,  SADAALSOR,  SORVERT,  and  HOLBOAT,  are  included  in  Appendix  C. 

It  is  reconmended  that  one  or  more  of  these  programs  be  modified  in  FY80. 

This  work  is  to  be  performed  primarily  by  NSWC,  (Mr.  J.  Renzi,  Dr.  P.  Huang, 
and  Mr.  J.  Matra). 

An  existing  computer  program  available  from  NSWC  Code  U14  will  provide  the 
hydroelastic  loading  for  the  SADAALSOR  program.  It  appears  that  the  SORVERT 
program  needs  no  such  mating  to  the  U14  program  since  it  already  contains 
hydroelastic  interaction.  However,  the  SORVERT  program  does  not  allow  oblique 
water  impact  where  the  SADAALSOR  program  has  no  such  restriction. 
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See  footnote  23  on  page  13. 
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Before  Che  dynamic  loading  case  can  be  considered,  a  shell  must  be  designed 
Co  withstand  hydrostatic  pressure  loading.  Such  a  design  constitutes  the  base¬ 
line  for  the  design  procedure  of  the  dynamic  case.  The  purpose  of  including  the 
HOLBOAT  program  is  to  have  a  good  closed  form  analysis  program  which  could  become 
a  design  program,  when  modified,  to  design  a  weight  efficient  shell  subject  to 
static  design  loads  such  as  hydrostatic  pressure.  The  HOLBOAT  program  is  the 
most  complete  and  versatile  program  available  in  terms  of  loading  and  wall 
construction.  It  is  also  in  closed  form  and  is  based  upon  the  proper  theories. 

A  complete  research  program  must  contain  some  method  of  verifying  the  accuracy 
of  the  analytical  model.  In  this  respect  the  NSWC  facilities  at  White  Oak  have 
the  capability  and  expertise  to  conduct  a  comprehensive  test  program  of  this  scope. 
The  Shock  Facilities  (E21)  conduct  water  entry  simulations  using  shock  tubes. 

The  Hydroballistics  Facility  (U14)  performs  scaled,  propelled  water  entry  tests. 

The  explosives  research  division  (R12)  has  the  capability  of  conducting  tests 
similar  to  those  of  Linberg  and  Herbert.^  Full  scale  model  tests  can  be  done 
at  the  NOSC  facility  at  Morris  Dam,  Los  Angeles.  The  type  and  extent  of  testing 
must  be  determined  by  considering  the  data  requirements  resulting  from  the 
parametric  study. 

Some  data  exists  which  may  be  used  in  its  present  condition.  In  FY77  a  water 
entry  program  was  conducted  in  which  a  lightweight,  high-strength  case  was  mated 
to  CAPTOR  flight  gear  and  air  dropped  twice  into  the  Charlie  Range  of  the  Patuxent 
River  at  the  NSWC  Solomon’s  facility.  The  case  was  of  sandwich  construction  with 
S-glass-epoxy  laminate  skins  and  an  aluminum  honeycomb  core.  Acceleration  and 
strain  data  were  gathered  and  film  coverage  was  provided.  This  data  will  provide 
one  excellent  means  of  comparison  with  analytical  results. 
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See  footnote  13  on  page  11. 
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Ref.  No.  1 

BRUSH,  D.  0.  and  MORTON,  F.  G.  May  1961 


"Stress  and  Stability  Analysis  of  Cylindrical  and  Conical  Shells." 


An  experimental  investigation  of  unstiffened,  thin-walled  cylindrical  shells 
subjected  to  transverse  loads  applied  by  a  cushioned  saddle,  this  report  does 
not  even  vaguely  resemble  the  kind  of  information  required  by  the  present 
investigation. 

Critique: 

1.  For  very  thin  shells. 

2.  For  surface  loads  induced  by  saddles. 

3.  Shell  treated  as  a  column  with  eccentric  load 

4.  Results  were  to  indicate  type  of  saddle  material  required. 

5.  All  results  based  on  experiment. 
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Ref.  No.  2 

LAKSHMIKANTHAM,  C.  and  GERARD,  G.  March  1964 


"Elastic  Stability  of  Cylindrical  Shells  Under  Axial  and  Lateral  Loads" 


Considers  the  forces  due  to  a  combination  of  axial  compression  and  external 
pressure  encountered  in  launch  vehicle  structures. 

"Using  the  Donnell  equation  for  small  deformations,  the  report  considers 
compressive  loading  combinations  on  the  stability  problem  of  an  unstiffened  circular 
cylinder."  (p.  773) 

The  Donnell  equation  in  terms  of  normal  displacement  in  the  absence  of  twisting 
forces : 


D  V^w 


+  V 


32w 

3x2 


+  N 


~2 
3  w 

3y2 


+  Et  9^w  =  0 

RZ  3x4 


Assumed  form  of  w:  w  =  sin  (x/Xx)  sin  (y/Xy)  where  Xx  and  Xy  are  the  half  wave- 
lengths  of  buckles  along  x  and  y  directions. 

Assume:  for  Ny  =  0,  Nx  =  N,  axial  load  only  Xx  =  L/tt  which  means  th^t  the  length 
of:  the  cylinder  *  *s  a  wavelength  of  deformation. 


Let : 

X  /X  =  B- 
x  y 

and 

T  -  Ny/N 

Then: 

N  ~  D  =  —  *  K-  +  EtL4/  (tt4R2D)  (1  +  02t)  (1  +  B2)2 

iT  (1  +  P  T ) 

Let  Z  =  (L2/Rt)  (1  -  V2)15 

kx  =  4(1  +  B2)2/(2  +  3B2t  +  t) 

where  is  the  buckling  stress  coefficient.  Then  plot  k-Z. 


Conclusions :  If  small  T  values  are  generated  as  a  result  of  initial  imperfections, 

or  prebuckling  deformations,  then  there  may  be  an  unusual  sensitivity  to  small 
values  of  T.  Buckling  stress  is  lowered  as  t  increases. 
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Ref.  No.  3 

SINGER,  BARUCH,  HARARI  May  1966 


"On  Stabili'.y  of  Eccentrically  Stiffened  Cylindrical  Shells  Under  Axial  Compression" 


This  paper  studied  the  effect  of  stiffener  eccentricity  on  the  critical  load 
for  cylindrical  shells  under  axial  compression.  The  behavior  of  the  eccentricity 
effect  depends  upon  the  geometry  of  the  shell  (height,  diameter,  and  thickness) 
while  the  geometry  of  the  stiffeners  only  influences  the  magnitude  of  this 
behavior. 

Only  2  of  the  8  possible  end  conditions  are  studied  assuming  that  secondary 
boundary  conditions  are  less  pronounced  in  stiffened  shells.  Therefore  only  simple 
supports  and  clamped  ends  are  examined. 

Assumptions : 


1.  The  stiffeners  are  distributed  over  whole  surface. 

2.  Normal  strains  vary  linearly  la  stiffener  and  sheet  and  are  equal  at  their 

point  of  contact. 

3.  Shear  membrane  force  carried  entirely  by  sheet. 

4.  Torsional  rigidity  of  stiffener  cross  section  is  added  to  that  of  sheet. 

Donnell  type  stability  equations  for  general  instability  initial  displacements 
are  assumed: 

1.  Similar  to  those  proposed  by  Batdorf  for  clamped  shell. 

2.  Which  solve  the  Donnell  type  equation  in  the  presence  of  boundary  conditions 

Rings  are  most  effective  in  stiffening  against  hydrostatic  pressure.  Stringers 
are  most  effective  in  stiffening  against  axial  compression.  Assumption:  The 
total  geometric  bending  stiffness  of  the  combined  stringer-shell  cross-section 
is  not  affected  by  the  position  of  the  stringers  (inside  or  out)  even  though  the 
moment  of  inertia  is  larger  for  outside  stringers.  For  closely  spaced  stringers 
this  difference  is  very  small  and  is  neglected.  Therefore,  by  "distributing" 
the  stringers  the  author  takes  care  of  both  inside  and  outside  types  because  they 
are  assumed  the  same.  _ 

Results:  The  two  modes  of  compressive  failure,  symmetric  and  antisymmetric  must 
be  considered  though  the  antisymmetric  mode  usually  yields  higher  buckling  loads 
except  for  very  long  and  thin  shells. 
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There  are  two  effects  of  eccentricity  of  stringers: 

1.  Primary  effect  -  outside  stringers  increase  the  extensional  stiffness 
in  the  longitudinal  direction  more  than  inside  stringers. 

2.  Secondary  effect  -  inside  stringers  increase  the  extensional  stiffness 
in  the  circumferential  direction  more  than  outside  stringers. 

When  subjected  to  axial  compression  rings  play  an  important  role  in  that  they 
act  to  shorten  the  effective  length  of  the  shell  and  thereby  put  the  deformation 
in  the  short  shell  regime.  "Now,  for  short  shells  the  main  resistance  to  buckling 
is  in  the  longitudinal  direction."  Experimentation  on  350  shells  was  conducted 
to  verify  these  calculations. 

Conclusions  included: 

1.  Outside  stringers  yield  higher  buckling  loads  than  inside  ones  for  all 
practical  geometries. 

2.  Eccentricity  effect  has  a  pronounced  maximum  at  values  which  are  common. 

3.  Rings  alone  are  much  less  efficient  as  stiffeners. 

4.  With  inside  rings  the  buckling  load  is  reduced  by  the  eccentricity,  whereas 
with  outside  rings,  the  axisymaetric  pattern  which  is  not  influenced  by  eccentri¬ 
city,  dominates. 

5.  Behavior  of  the  eccentricity  effect  depends  on  the  geometry  of  the  shell 
while  the  stiffeners  influence  its  magnitude. 

Critique: 

1.  For  static  case  only,  otherwise,  much  data  on  stiffeners  giving  much  insight 
into  failure  modes  and  stiffener  influence. 

2.  Used  small  deflection  theory. 
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HEGEMIER,  G.  A.  January  1967 


"Stability  of  Cylindrical  Shells  Under  Moving  Loads  By  the  Direct  Method  of 
Liapunov" 


Investigates  the  stability  of  a  long,  thin,  elastic,  circular  cylindrical 
shell  subjected  to  axial  compression  and  an  axisymmetric  load  moving  with  constant 
velocity  along  the  shell  axis. 

Governin£jS<juation£:  Direct  method  of  Liapunov  using  nonlinear  Donnell-type  shell 
theory. 

Critique: 

1.  Analytic  only 

2.  Non  dynamic  analyses 

3.  No  design  criteria 

4.  Steady  state  response  only 
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Ref.  No.  5 

HUTCHINSON,  J.  W.  and  AMAZIGO,  J.  C.  March  1967 


"Imperfection  Sensitivity  of  Eccentrically  Stiffened  Cylindrical  Shells" 

This  paper  is  a  quantitative  study  of  the  effect  of  stringer  and  ring  stiffening 
on  the  sensitivity  of  circular  cylindrical  shells  to  imperfection  while  subjected 
to  axial  compression  and  hydrostatic  pressure. 

Governing  Equations:  Nonlinear  Donnell-type  strain  displacement. 

Assumgtion:  Imperfection  is  assumed  in  the  shape  of  the  buckling  mode. 

The  paper  goes  into  detailed  analysis  of  qualitative  effects  of  rings 

2  k 

and  stringers  on  imperfection  sensitivity  with  respect  to  Z,  Z  *  ^(1-v  )  . 

Results:  Under  axial  load,  for  axially  stiffened  shells,  the  imperfection  sensiti¬ 
vity  and  buckling  load  are  strongly  dependent  on  whether  stringers  are  inside 
or  outside  the  cylinder.  Stiffening  can  reduce  or  eliminate  imperfection 
sensitivity. 

Critique: 

1.  For  static  buckling  case  only. 

2.  No  design  information,  not  general. 

3.  Nonlinear  theory. 
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Ref.  No.  6 

HORTON  and  TENERELLI  December  1967 


"The  Instability  of  Unstiffened  and  Ring  Stabilized  Thin  Walled  Circular  Cylinders 
Under  Non-Uniform  Axial  Load  Conditions" 


This  work  deals  with  the  aero-space  industry  problems  and  therefore  includes 
temperature  effects.  After  a  rather  complete  review  of  the  previous  works  on 
axial  and  temperature  induced  loading  of  right  circular  cylindrical  shells,  tests 
are  performed  to  eliminate  inconsistency  in  the  findings. 

Their  findings:  analysis  shows  that  critical  stress  should  be  independent 
of  distribution  whereas  experimentation  shows  this  is  not 
true. 

The  critical  pure  compressive  stress  is  *  30-4031  less  than  the  critical  pure 
bending  stress  experimentally. 

Two  experiments: 

1.  Beverage  cans. 

2.  Elastic  buckling  using  interior  mandrel  to  restrict  buckling  amplitude. 

Analysis  was  accomplished  using  a  computer  program  which  solved  the  three  simul¬ 
taneous  equations  of  Fliigge  (displacement  equations). 

"The  paper  shows  that  both  perfect  and  imperfect  unstiffened  thin- 
walled  right  circular  cylindrical  shells  buckle  under  non-uniform  load 
conditions  when  the  maximum  compressive  stress  reaches  the  value 
which  would  cause  instability  under  uniform  load  conditions."  (p.  23) 

Therefore,  for  the  uniform  imperfection,  the  design  stress  level  should  be  the 
same  for  the  uniform  compressed  states  as  for  the  varying  axial  stress  condition. 

But  due  to  imperfection  nonuniformity,  increased  allowable  stress  may  be  used 
in  flexure  provided  it  is  possible  to  assess  the  probability  that  the  location 
of  the  weakest  section  of  the  shell  will  not  correspond  to  the  region  of  maximum 
stress. 

Critique: 

1.  Unstiffened  shells  only. 

2.  Nonlinear  analysis. 

3.  Imperfections  dealt  with  as  an  aside. 
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Ref.  No.  7 

SHUMIK,  M.  A.  (Moscow)  June  1974 


"The  Stability  of  Cylindrical  Shells  Under  the  Action  of  Dynamic  Radial  Pressure" 

The  study  determines  that  there  is  a  critical  pressure  impulse  for  a  given 
shell . 

Assumgt^ions : 

orthotropic 
ring  stiffened 

disregard  shear  &  torsional  rigidity  of  ribs 

loading  is  not  too  "fast",  therefore  can  disregard  the  tangetial  forces 
of  inertia  and  the  subcritical  stressed  state  is  inertialess 

initial  and  subsequent  deflections  assumed  sinusoidal. 

A  Conclusion: 


Critical  pressure  does  not  rely  upon  the  cylinder's  length. 

Critique : 

1.  Very  hard  to  read. 

2.  Radial  pressure  dynamic  load. 

3.  Very  slow  pressure  pulse  is  assumed  so  that  inertia  may  be  neglected  which 
makes  this  a  "quasi-static"  problem,  not  a  dynamic  problem. 
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Ref.  No.  8 


EL  NASCHIE ,  M.  S. 


June  1975 


"Localized  Diamond  Shaped  Buckling  Patterns  of  Axially  Compressed  Cylindrical 
Shells" 


This  paper  proposes  a  mathematical  method  of  analyzing  a  shell  based  on  a 
two-dimensional  Riemann  space  model  interpreted  as  the  buckling  of  a  strut  on 
a  "bending"  foundation  with  an  attenuated  eigenvector.  The  method  uses  a  simple 
differential  equation  with  Pogorelov  coefficients  to  result  in  a  value  of  the 
lower  critical  stress  (oc)  for  a  given  buckling  mode  (w) 

0  =  0.093  - — — 

C  r(l-U2) 


2ae  £(3)**  cos  6x  -  sin  6xJ 


where  y  and  6  are  functions  of  thickness  (t),  radius  (r)  and  the  second  fundamental 
tensor  (b): 


y  =  K/(2)li, 


K 


(30tVr2) 

(bW)) 


* 


Critique : 

The  schema  is  different  but  does  not  pursue  the  subject  in  the  desired  direction. 
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Ref.  No.  9 


WELLER,  T.  and  SINGER,  J. 


December  1977 


"Experimental  Studies  On  the  Buckling  Under  Axial  Compression  of  Integrally  Stringer- 
Stiffened  Circular  Cylindrical  Shells" 


This  paper  is  a  summary  of  other  works  by  the  authors. 

Panels  were  designed  using  criterion  derived  by  Koiter  in  1956. 

General  instability  calculated  with  "smeared"  stiffener  linear  theory. 
Conclusions : 

1.  Linear  theory  can  be  used  effectively  as  a  first  approximation  for  predic¬ 
ting  buckling  loads  of  stringer  stiffened  shells. 

2.  Stringer  stiffened  shells  are  sensitive  to  boundary  conditions. 

3.  Agreement  between  experiment  and  theory  relies  primarily  on  stiffener 
spacing,  stiffenei  area  ratio,  and  shell  geometry  parameter  Z. 

4.  Material  nonlinearity  may  have  significant  affect  on  linear  theory 
applicability. 

5.  Stringer  stiffened  shells,  within  the  range  of  geometry  in  practical  design, 
are  structurally  more  efficient  than  equivalent  weight  isotropic  shells. 

Critique :  Static  Case  -  eccentricity  of  loading  avoided.  Linear  theory  was 
compared  with  the  experimentation. 
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Ref.  No.  10 

YAO,  J.  C.  December  1961 

"Stability  of  a  Cylinder  Under  Dynamic  Radial  Pressure" 

This  paper  investigates  the  stability  of  an  infinite  cylinder  under  a  uniform 
radial  pressure  applied  as  a  rectangular  shaped  function  of  time. 

Motion  is  assumed  to  consist  of  a  uniform  radial  displacement  plus  one  flexural 
mode  having  two  waves  around  the  circumference. 

Coupling  -  axisymmetric  mode  with  flexural  mode. 

Results :  Flexural  mode  was  shown  to  respond  in  an  oscillatory  or  exponential 
manner  depending  on  the  magnitude  of  the  pressure. 

Critique : 

1.  End  conditions  assumed  trivial  to  case  examined. 

2.  No  stiffening  of  any  kind. 

3.  Not  general. 

4.  Assumptions  made  to  ignore  inertia  terms. 
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Ref.  No.  11 

YAO,  J.  C.  June  1963 

"Dynamic  Stability  of  Cylindrical  Shells  Under  Static  and  Periodic  Axial  and 
Radial  Loads" 

Assumptions : 

1.  The  load:  one  static  and  one  periodic  load  component  only. 

2.  The  frequency  of  the  periodic  load  is  «  the  frequency  of  stress  waves 
in  the  cylinder. 

3.  Radial  displacement  of  the  form  w  =  £(t)  sin  not  cos  k$ . 

4.  Loading  expression  assumed  to  be  Nj  =  Nq  (a  cos  wT  +  e) 

N2  =  rP  =  rPQ  (b  cos  wT  +  f) 

therefore  defining  periodic  end  loads  and  radial  pressures. 

5.  Loads  have  identical  frequencies. 

6.  Both  ends  are  simply  supported. 

Using  Mathieu's  equation  and  a  modified  Donnell's  equation  of  motion  stability 
of  the  system  with  the  above  assumptions  is  examined  by  seeking  two  wave  parameters. 

Individual  cases  are  examined  by  making  other  assumptions.  Of  the  cases 
examined,  one  includes:  Cylinders  Under  Static  &  Periodic  Axial  Loads. 

Assumed :  Cylinders  are  assumed  to  be  of  medium  length  so  end  conditions  have 
no  influence  on  their  behavior. 

Results :  When  w/wq  becomes  larger,  stability  region  increases. 

Critique : 

1.  End  conditions  ignored,  therefore  not  very  general. 

2.  No  stiffening  of  any  kind. 

3  Sinusoidal  loading  only. 

4.  Form  of  displacement  assumed. 
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Ref.  No.  12 

ROTH  and  KLOSNER  October  1964 


"Nonlinear  Response  of  Cylindrical  Shells  Subjected  to  Dynamic  Axial  Loads" 


This  paper  investigates  the  stability  of  circular  cylindrical  shells  subject 
to  uniformly  distributed,  suddenly  applied  axial  loads  of  constant  magnitude  of 
prolonged  or  finite  time  duration.  Analysis  includes  a  4  degree  of  freedom  system 
studied  numerically. 


Assumptions : 

1.  Longitudinal  inertia  neglected. 

2.  Effect  of  edge  restraints  negligible. 

3.  Loading  in  the  form  of  the  heavyside  step  function. 

Method :  Use  Donnell's  strain-displacement  relationships  and  Hooke's  law  for  thin 
shells  and  an  equation  for  elastic,  extensional  strain  energy.  Then  the  equations 
of  motion  are  found  by  applying  Hamilton's  principle.  Now  introduce  an  Airy  stress 
function.  Write  equilibrium  and  compatibility  equations  to  give  2  nonlinear 
differential  equations  to  be  solved  for  stress  function  F  and  radial  displacement 
w. 

To  solve  these: 

-  approximate  radial  deflection  function 


w  =  h  jaj(t)  cos  ^  cos  ^  +  a^ft)  cos  +  a^Ct)  cos  — +  a^(t) 
L  X  y  x  y 


-  initial  imperfections  form 


—  ,  ,,  ttx  Try  .  ,  2hxn 

w  -  h(dj  cos  cos  +  “2  cos  TT"' 
x  y  x 


Then  substitute  into  equilibrium  and  compatibility  equation  and  apply  Hamiltons 
principle.  This  yields  4  simultaneous  nonlinear  differential  equations. 


These  are  solved  with  8  initial  conditions  for  5  parameters  using  a  Runge- 
Kutta  formula. 


Results :  A  series  of  plots  showing  solutions  to  the  radial  displacement  function 
versus  time  for  given  load  conditions.  For  zero  imperfections  the  dynamic 
buckling  stress  coincided  with  that  of  the  classical  linear  theory.  The  most 
important  factors  influencing  peak  resistance  are  the  initial  imperfections, 
initial  stresses,  and  deviations  from  uniform,  isotropic-elastic  behavior.  When 
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a  small  initial  displacement  or  a  small  initial  velocity  was  assumed,  a  sudden 
jump  in  the  amplitude  of  the  ai  mode  did  not  occur  until  the  stress  reached  a 
value  equal  to  the  classical  buckling  stress.  Initial  imperfections  play  a  dominant 
role  in  the  reduction  of  peak  stresses  achieved  prior  to  buckling.  The  dynamic 
buckling  stress  a^so  decreases  with  increase  in  time  duration  of  loading.  The 
critical  impulse  Ig  =  Og  X  decreases  as  the  time  duration  of  loading  decreases. 

Conclusions :  Many  and  varied,  too  many  to  mention  here. 


Critique : 

1.  All  inertia  terms  neglected,  therefore  raising  doubts  as  to  accuracy  of 
model  and  applicability  to  our  problem. 

2.  By  assuming  the  modes  of  deflection  and  imperfection  the  model  is  not 
as  general  as  may  be  desired. 

3.  The  model  was  for  right  circular  cylinders  and  therefore  not  of  the  general 
nature  most  desirable. 

4.  It  did  however  help  define  governing  parameters  and  give  an  idea  of  the 
directions  to  pursue. 
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Ref.  No.  13 

LINDBERG,  H.  E.  and  HERBERT,  R.  E.  March  1966 


"Dynamic  Buckling  of  a  Thin  Cylindrical  Shell  Under  Axial  Impact" 


This  paper  examines  buckling  of  thin  cylindrical  shells  from  axial  impact 
assuming: 

1.  Initial  imperfections  can  be  approximated  by  "white  noise". 

2.  Linear  small  deflection  theory  applies. 

3.  Simplified  Donnell-type  equations  give  satisfactory  accuracy. 

4.  Shell  is  simply  supported. 

5.  Imperfections  and  final  buckled  form  are  stationary  (in  space). 

Theory  is  used  to  calculate  the  resulting  growth  of  normal  displacement  modes. 
By  using  linear  theory  the  buckling  form  is  allowed  to  assume  its  "own"  shape 
while  deflections  are  small.  A  statistical  analysis  gives  the  expected  values 
for  the  "preferred"  axial  and  circumferential  wavelengths. 

Experimentation:  Run  simultaneously  with  theoretical  investigation  to  see  if 
the  shell,  once  deformed  into  its  linear  shell  theory  initial  deflection,  would 
continue  in  that  form.  Methods  included: 

1.  Shell  free  at  end  opposite  impact. 

2.  Impacted  end  bonded  to  massive  inside  ring  to  provide  a  clamped  boundary 
to  the  shell. 

3.  End  ring  accelerated  explosively  so  that  time  and  simultaneity  of  impact 
could  be  controlled  to  within  ~  2  microseconds. 

4.  High  speed  cameras  used  240,000  frames/sec. 

Conclusions :  Experimentation  confirmed  the  premise  that  early  small-deflection 
buckling  dictates  the  pattern  into  which  large-deflection  buckles  form. 

For  calculating  the  threshold  of  buckling  from  axial  impact,  it  may  be  suffi¬ 
cient  merely  to  define  some  magnitude  of  the  linear-amplification  function  as 
indicative  of  buckling. 
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Critique : 

1.  Uniform  axial  load  is  not  the  same  as  our  problem  in  general,  but  close. 

2.  Linear  theory. 

3.  Initial  imperfections  included. 

4.  Simply  supported  shell,  therefore  assuming  only  one  kind  of  end  condition. 

Later  paper  shows  that  for  high  speed  impact  (impact  >  V  )  the  buckling  was 
axisymmetric .  c 
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Ref.  No.  14 

COPPA,  A.  P.  July  1966 

"EffecC  of  End  Conditions  on  the  Buckling  of  Cylindrical  Shells  Under  Axial 
Compression  Impact" 


This  is  a  report  on  experimentally  observed  effects  of  edge  conditions  on 
the  buckling  behavior  of  cylindrical  shells  under  axial  impact. 

Experiments  were  performed  using  cylinders  aet  on  a  stand,  pressurised  with 
a  free  end  and  impacted  by  a  drop  test  apparatus. 

Tests  include  different  pressures  and  end  restraints. 

Results : 

1.  Increasing  pressure  increases  buckling  stress  (to  .799....). 

2.  Lateral  restraint  increases  buckling  stress  (to  .860  classical  static). 

3.  Without  restraint  against  inward  displacements,  local  collapse  at  the 
impacted  end  occurs  more  rapidly. 

4.  Asymmetric  impact  precludes  initial  buckling  response  characteristic  of 
symmetrical  impact.  Only  small  stress  waves  were  transmitted  to  midlength  posi¬ 
tions  and  local  buckling  occurred  at  the  position  first  struck  by  the  impact  head. 

5.  Buckling  stress  increases  with  velocity. 

Insights :  Triangular  buckles  represent  a  lower  energy  mode  of  response  than  the 
symmetrical  mode.  Symmetrical  buckles  result  from  predominantly  outward  displace¬ 
ments,  (triangular  buckles  inward).  When  a  shell  is  impacted,  it  is  forced  to 
respond  quickly  in  an  energetically  economic  deformation  mode.  Therefore  the 
dominant  mode  is  axisymmetric  initially.  Whether  this  mode  prevails  depends 
upon:  shell  thickness,  amount  of  internal  lateral  support,  relationship  between 
impacted  strain  and  the  yield  strain  of  the  material. 

Conclusion:  The  triangular  response  mode  is  the  mode  preferred  away  from  the 
end  restraint,  and  is  independent  of  edge  conditions. 


Critique : 


1.  No  mathematical  analysis  (theory). 


2.  No  general  quantitative  design  criteria. 


3.  Qualitative  only. 
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Ref.  No.  15 


DIETZ,  W.  K. 


January  1967 


"On  The  Dynamic  Stability  of  Eccentrically  Reinforced  Circular  Cylindrical  Shells" 


An  analytical  investigation  using  large  deflection  theory  including  initial 
imperfections  and  rotatory  inertia  dealing  with  dynamic  buckling  loads  of  a  clamped 
reinforced  circular  cylindrical  shell. 

Assumptions: 

1.  Stringer  and  ring  stiffeners  are  closely  spaced  so  that  an  "equivalent 
shell"  approach  may  be  used  where  the  reinforcement  effects  are  smeared  out  over 
the  whole  surface. 

2.  Radial  displacement  is  assumed  to  be  diamond  patterned  and  satisfy  clamped 
boundary  conditions. 

-yt 

3.  End  shortening  of  the  form  V  *  VQ  e  . 

Conclusions: 

1.  Rotatory  inertia  can  be  neglected  in  determing  the  critical  dynamic  buckling 
load. 

2.  The  size  of  the  imperfection  amplitude  affects  the  critical  dynamic  buckling 
load  drastically.  Whether  reinforcement  decreases  this  sensitivity  or  not  is 

not  addressed. 

3.  The  effect  of  imperfection  direction  is  unimportant. 

4.  The  exponentially  decaying  end  shortening  reduces  dynamic  buckling  load 
which  is  dependent  on  the  time  constant  1/y. 

Critique: 

1.  Done  for  one  particular  shell,  not  general  design  data. 

2.  Dynamic  case. 

3.  Large  Deflection  Theory. 
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Ref.  No.  16 

CROMER,  C.  C.  April  1969 


"An  Investigation  of  the  Nonlinear  Dynamic  Response  of  Cylindrical  Shells  Under 
Transient  Pressure" 


A  ring  stiffened  nearly  circular  cylindrical  shell  of  finite  length  under 
transient,  axisymmetric  radial  loads  of  arbitrary  axial  distribution. 

Governing  Equation:  Nonlinear  Donnell-type  solved  using  Fourier  series  expansions 
of  the  dependent  variables.  Newmark's  beta-method  combined  with  Gauss  elimination 
for  time  integration. 

This  paper  was  written  for  R.  E.  Ball,  to  model  the  Dahlgren  shock  tube  under 
operational  loading  conditions. 

Critique: 

1.  Ring  stiffened  only. 

2.  Axisymmetric  radial  loads  only. 

3.  Inertia  not  considered. 

4.  Initial  imperfections  symmetric  about  a  plane  thru  the  axis  are  accounted 

for. 

5.  Nonlinear  coupling  between  modes  retained. 

6.  Boundary  conditions  consistently  observed. 

7.  Effects  of  discrete  eccentric  rings  included. 

8.  Viscous  damping  included. 
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Ref.  No.  17 


MORTIMER,  R.  W. ,  CHOU,  P.  C.  and  ROSE,  J.  L. 


October  1969 


Annual  Report  for  Wave  Propagation  In  Stepped  and  Joined  Shells 


A  computer  program  (MCDIT-21)  was  utilized  to  analyze  shell  impact  problems. 
Impacts  were  characterized  as  longitudinal  velocity  step  functions. 

Tests  were  conducted  using  pendulum  impactor,  photodiodes  and  light  source. 
They  verified  the  program's  output. 

A  modified  membrane  theory  is  used.  If  differs  from  the  classical  theory 
in  that  transverse  shear  deformation  is  included. 

Critique; 


1.  No  inertia  effects. 

2.  Simple  unstiffened  shells. 

3.  No  imperfection  effects  included. 
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Ref.  No.  22 


DAREVSKII,  V.  M. 


"Stability  of  A  Cylindrical  Shell  to  Dynamic  Axial  Loading" 

This  article  examines  the  loss  of  stability  by  a  cylindrical  shell  for 
several  types  of  dynamic  loading.  The  case  of  instantaneous  axial  compression 
is  one  of  these. 

Instability  is  defined  as  the  displacements  of  points  in  the  elastic  body 
with  no  change  in  their  velocities  or  accelerations.  This  occurs  at  a  time 

c  =  V 

The  basic  equations  used  are  the  nonlinear  equations  of  motion  of  the 
elastic  body  with  allowance  for  initial  deformation: 

A^ (o^  i  O 2 >  •  •  •  )  ^  Oj i  •  •  •  5  v,  w  )  ^  t  q^ »  q^  5  w) 


^  32Ui 

+  q,  -  $  — 7  ■  0 

3t^ 


<i  *  1,  2,  3) 


where  6  =  material  density 

t  =  time 

U1  *  u>  U2  =  v»  U3  *  w 
q^  =  loads 

One  case  investigated:  medium  length  cylindrical  shell 

one  edge  clamped,  the  other  loaded 

the  load  is  evenly  applied  over  the  entire  end. 

The  case  of  instantaneous  axial  compression  is  examined  using  a  shell  of 
"medium"  length,  with  one  edge  clamped  and  the  other  under  uniform  axial  compressive 
force  -P,  instantaneously  applied. 

The  critical  load  P  is  determined  using  simplified  secondary  displacement 
equations  analogous  to  the  stability  equations  for  static  forces. 

The  aim  is  to  find  the  relative  value  of  Pc  to  the  corresponding  critical 
static  load  (in  its  classical  sense).  After  much  analysis: 

1  pc8  r^~k  -  pc  <  i  pc8  min  (i»  b 


where: 


3v  (1  -  v  (1  +  cos  2p)  _1  nr 
2V2W  R*h 
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No.  18 

BALL,  R.  E.  1972 


"A  Program  for  Che  Nonlinear  SCaCic  and  Dynamic  Analysis  of  Arbitrarily  Loaded 
Shells  of  Revolution" 


This  paper  discusses  a  digital  computer  program  which  analyzes  the  response  of 
arbitrarily  loaded  shells  of  revolution. 

Governing  Equation;  Sander's  non-linear  thin  shell  theory  for  conditions  of  small 
strains  and  moderately  small  rotations. 

Dependent  variables  are  approximated  using  Fourier  series.  Nonlinear  modal 
coupling  is  treated  with  pseudo  loads.  Derivatives  approximated  using  central 
finite  differences.  Displacement  accelerations  by  implicit  Houbolt  backward 
difference. 

Program  used  ONLY  when: 

1.  Geometries  and  material  properties  are  axisymmetric. 

2.  Applied  pressure  and  temperature  distributions  can  be  expressed  in 
a  Fourier  sine  or  cosine  series. 

3.  Isotropic  material,  V  *  constant,  E  may  vary  through  thickness. 

4.  Boundaries  may  be  closed,  free,  fixed  or  restrained. 

Critique: 

1.  Cannot  treat  stringers  or  discrite  rings. 

2.  Demands  interpretation  of  results  to  ascertain  whether  shell  has  reached 
its  buckling  load. 

3.  Does  not  include  initial  imperfections. 

4.  Nonlinear  theory. 

5.  Versatile  boundary  conditions. 
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R  and  h  are  the  radius  and  thickness  of  the  shell,  and  E  and  V  are  the  elastic 
modulus  and  Poisson  ratio,  respectively. 

Critique: 

1.  For  "short”  shells  with  thickness  small  so  they  can  be  classified 
as  "medium"  length. 

2.  No  experimentation. 

3.  No  stiffness  accounted  for. 

4.  No  imperfection. 

5.  Boundary  conditions  satisfied. 

6.  Nonlinear  theory. 
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Ref.  No.  23 


BHATIA,  P.  and  BABCOCK,  C.  D.  Jr. 


February  1973  I 

i 


"Dynamic  Buckling  of  Structures" 


This  report  uses  an  energy  approach  to  examine  the  stability  of  arches 
and  spherical  caps  subject  to  dynamic  loads. 

Definition:  The  Dynamic  Buckling  Load  is  defined  as  the  scalar  measure  of  that 
loading  history  at  which  there  is  a  finite  jump  in  an  appropriate  measure  of  the 
response  for  infinitesimal  increase  in  the  load  measure. 

Critique:  Only  two  degree  of  freedom  systems  are  examined  though  it  is  mentioned 
that  the  theory  can  be  applied  to  higher  degree  systems.  There  is  one  mention 
of  an  applicability  problem,  the  theory  must  not  apply  to  all  structures. 

Also,  complex  structures,  imperfections,  and  real  material  properties  pose 
problems  to  successful  use  of  this  approach. 
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Ref.  No.  24 


TAMURA,  Y .  S. 


June  1973 


"Dynamic  Stability  of  Cylindrical  Shells  Under  Step  Leading" 


This  paper  studies  dynamic  stability  of  imperfect  circular  cylindrical  shells 
under  axial  step  loading.  The  study  includes  effects  of  wave  number  of  radial 
mode  shape,  mass  on  the  loaded  edge,  and  damping  of  the  axial  motion. 

-  Axial  inertia  approximated 

-  Radial  inertia  included 

-  Radial  displacement  approximated  by  a  finite  degree  of  freedom  system 

-  Initial  imperfections  given  by  experimental  results 

-  Structural  damping 

-  Critical  loads  determined  by  numerical  integration  of  eq.  of  motion 

-  Nonlinear  Donnell's  equations. 


Method; 

1.  Assume  displacement  function  w  and  imperfection  function  w. 

2.  Determine  inplane  displacements  from  their  equation. 

3.  Evaluate  total  Potential  &  Kinetic  energy. 

4.  Apply  Hamilton’s  principle. 

Results:  Gives  proper  imperfection  sensitivity  a  4  degree  of  freedom  model. 

Conclusions :  Through  parametric  studies  of  the  wave  number,  mass,  and  damping 
factor,  due  to  frequency  coupling  between  axial  and  radial  motions,  the  axial 
inertia  plays  an  essential  role  in  characterizing  the  dynamic  instability  of  a 
finite  length  shell. 

Damping  suppresses  parametric  resonance  between  radial  and  axial  modes  and 
can  bring  the  buckling  load  close  to  the  static  value. 

Critique: 

1.  No  stringer,  ring  inclusion. 

2.  Inertia  terms  included  (axial  and  radial). 

3.  Damping  included. 

4.  Nonlinear  theory. 

5.  Initial  imperfections  included. 
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Ref.  No.  25 


LANGE,  et.  al. 


1973 


"The  Role  of  Resonant  Interactions  in  Dynamic  Stability  of  Elastic  Systems" 


This  is  a  mathematical  treatise  dealing  with  the  resonance  phenomenon  in  elastic 
systems.  The  dynamic  stability  of  a  cylindrical  shell  under  axial  compression 
illustrates  the  work.  The  load  can  be  step  or  impulsive. 

Basic  Equations:  Nonlinear  von  Karman-Donnell 

The  analysis  includes  radial  inertia  force.  Work  is  done  directly  with  the 
equations  for  radial  displacement  w  and  the  Airy  stress  function  f  representing 
deviations  from  the  membrane  prebuckling  solution  for  the  perfect  circular 
cylindrical  shell.  Time  is  made  nondimensional : 


(same  as  Tamura  and  Bobcock,  March  1975). 

* 

where  t  is  real  time. 

Results :  Nonlinear  coupling  of  normal  modes  of  the  solutions  to  the  equilibrium 

equations  can  have  significant  bearing  on  the  stability  of  elastic  systems.  In 
general  the  analysis  investigates  the  phenomenon  of  resonance  and  its  importance 
to  system  stability.  It  does  not  delve  any  further  into  the  subject  of  shells 
than  already  stated. 

Critique : 

1.  No  stiffeners. 

2.  No  imperfection  analysis. 

3.  Theoretical  only,  no  comparison  with  experimentation 

4.  Nonlinear  theory. 
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Ref.  No.  26 

FISHER,  C.  A.  and  BERT,  C.  W.  November  1973 

"Design  of  Crashworthy  Aircraft  Cabins  Based  on  Dynamic  Buckling" 

Also:  Fisher,  C.  A.  and  Bert,  C.  W. ,  "Dynamic  Buckling  of  an  Axially  Compressed 
Shell  with  Discrete  Rings  and  Stiffeners",  Journal  of  Applied  Mechanics, 
1973. 


The  authors  investigate  the  structural  integrity  of  aircraft  structures  using 
nonlinear  analysis  developed  to  treat  dynamic  buckling  of  a  circular  cylindrical 
shell  having  discrete  rings  and  stringers  subjected  to  rapidly  applied  axial 
compressive  load.  Both  general  or  localized  instability  are  considered. 

Critique : 


1.  Program  uses  tremendous  amount  of  computer  time  and  core  space. 

2.  Outside  stiffeners  treated  perfunctorily. 

3.  Imperfections  neglected. 
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Ref.  No.  27 

TSUI  and  LAKSHMIKANTHAM  June  1973 


"The  Nonlinear  Dynamic  Response  of  Stiffened  Shells  Under  Compressive  Axial  Impulse"  j 

i 


The  theoretical  model  includes: 

-  imperfections 

-  eccentricity  effects 

Equations :  Coupled  set  of  4th  order  nonlinear  von  Karman-Donnel 1  type.  Integrated 
using  4th  order  Runge-Kutta  scheme. 

Assumptions : 

1.  5-parameter  modal  representation  of  buckling  pattern. 

2.  Asymmetric  initial  imperfections. 

Results : 


1.  Shows  that  heavily  stiffened  shells  are  not  insensitive  to  imperfections 
in  the  dynamic  loading  case. 

2.  Dynamic  buckling  loads  are  less  with  inside  stiffeners  (similar  to  static 
studies)  such  as  Hutchinson  and  Amazigo) . 

3.  With  an  increase  in  imperfections,  the  effects  of  stiffener  eccentricity 

may  vanish,  (increasing  imperfection  size  decreases  the  differences  made  by  stiffener 
eccentricity  (+  or  -)  to  where  eccentricity  effects  may  vanish  for  the  axial  case 
and  decrease  considerably  for  the  general  case. 

Critique:  Nonlinear  theory  -  The  parametric  analysis  results  from  an  involved 
computer  program  which  numerically  integrates  four  nonlinear  differential  equations 
(time  steps  using  tremendous  computer  time).  No  inertia  terms.  These  people 
are  on  the  right  track. 
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Ref.  No.  28 


TSUI,  et.  al. 


December  1973 


"Linear  Stability  of  Eccentrically  Stiffened  Cylindrical  Shells  Under  Axial 
Compression" 


This  paper  considers  the  effect  of  stiffener  eccentricity  on  buckling  modes 
rather  than  buckling  loads.  The  information  is  significant  in  that  the  knowledge 
of  buckling  modes  allows  the  theoretician  to  have  a  priori  knowledge  of  the  shape 
of  the  shell  failure. 

The  mathematics  deal  with: 

-  Stringer  stiffeners,  inside  and  out 

-  Ring  stiffeners,  inside  and  out 

Theory  Used:  Linear  theory  of  Singer  which  accounts  for  coupling  between  bending 
stress  resultants  and  middle  surface  stretching,  and  membrane  stress  resultants 
and  changes  in  the  middle  surface  curvature. 

-  Stiffeners  are  smeared 

-  Stiffeners  are  essentially  beam  like  elements 

-  Local  instabilities  ruled  out. 

Constitutive  equations  given  in  matrix  form. 

Critique:  This  kind  of  thing  may  be  useful  in  the  assumption  of  a  preferred  buckling 

mode  for  a  given  shell  type. 
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Ref.  No.  29 


LAKSHMIKANTHAM  and  TSUI 


February  1974 


"Dynamic  Stability  of  Axially  Stiffened  Imperfect  Cylindrical  Shells  Under  Axial 
Step  Loading" 


Investigates  the  dynamic  general  instability  of  stiffened  shells  under  a 
variety  of  impulsive  loads  and  deals  with  the  case  of  an  axially-stiffened 
(stringer-stiffened)  cylinder  under  axial  step  loading. 

Donnell-type  formulation  valid  for  moderate  length  cylinders 

2  4 

Z  =  ~  10  <  Z  <  10 


Derivation  is  similar  to  that  of  Baruch  and  Singer  except  for  inclusion  of 
strain-displacement  relationships  and  initial  imperfections. 

-  Stiffeners  are  close  enough  together  to  be  smeared  out  to  calculate 
effective  rigidities. 

-  Beamlike  elements  carrying  no  shear. 

-  Local  instabilities  are  ruled  out. 

Nonlinear  strain-displacement  relationship  and  initial  imperfections  are 
included. 

Conclusions;  For  dynamic  buckling: 

1.  Outside  stiffeners  are  imperfection  sensitive. 

2.  Imperfection  sensitive  structures  under  static  load  will  display  the  same 
property  under  step  loads. 

Equation  of  motion  -  same  as  for  Lakshmikantham  and  Tsui,  1975. 

Initial  Imperfections:  a  four-parameter  modal  approximation: 

w  »  (— )  (dj  cosfm  it  x)  cos(ny)  +  62  cos(2m  IT  x)) 


IT  =  IT  R/L 
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Ref.  No.  30 


LAKSHMIKANTHAM  and  TSUI 


September  1975 


"Dynamic  Buckling  of  Ring  Stiffened  Cylindrical  Shells" 


This  paper  is  in  the  same  vein  as  previous  papers  by  the  authors. 

1.  Donnell 's  nonlinear  shell  theory  is  applicable. 

?.  Stiffeners  are  close  together  and  therefore  effective  rigidities  are 
calculated  by  smearing. 

3.  Stiffeners  are  beam-like  elements  carrying  no  shear. 

4.  No  local  instability  failures  -  entire  shell  including  stiffeners 
is  activated  at  buckling. 

Longitudinal  and  tangential  inertial  terms  are  of  lower  order  of  importance 
compared  to  normal  inertia. 

Buckling  Criteria:  A  jump  in  the  peak  amplitude  of  the  deflection  is  associated 
with  a  critical  load. 

Equations  of  Motion;  Donne 11-Karman  type  coupled  PDE  derived  from  the  Hamilton 
Principle  followed  by  Galerkin-type  procedure.  Use  a  4  parameter  modal  approxima¬ 
tion,  representing  diamond-shaped  buckles. 

The  paper  investigates  axial  and  lateral  step  loads. 

Results:  Ring  stiffened  shells  are  imperfection  sensitive.  The  essential  mechanism 
for  imperfection  sensitivity  in  the  dynamic  range  is  the  so-called  "column" 
asymmetric  mode  with  m  =  1  and  n  =  1. 

Conclusions: 

1.  Axially  and  ring  stiffened  cylinders  are  imperfection  sensitive. 

2.  Internal  ring  stiffened  shells  are  imperfection  insensitive  under  axial 
step  loading.  The  addition  of  small  amounts  of  lateral  pressure  reverses  this 
effect. 


3.  "Ring  stiffened  shells  under  pure  lateral  pressure  are  imperfection 
sensitive  for  all  ranges  of  ring  heaviness." 

4.  In  the  dynamic  range  the  mechanism  of  imperfection  sensitivity  is  the 
column  asymmetric  mode. 

5.  Internal  rings  force  the  shell  response  into  wave  patterns  conforming 
to  ring  spacing  under  axial  load.  Under  a  lateral  load  this  effect  is  weakened 
and  eventually  produces  the  column  failure  mode. 
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Critique: 

1.  Stiffeners  smeared. 

2.  Longitudinal  and  tangential  inertia  deamed  unimportant 

3.  Nonlinear  theory  used. 

4.  Step  loads  are  examined.  This  may  not  be  realistic. 
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Ref.  No.  31 


TAMURA  and  BABCOCK 


March  1975 


"Dynamic  Stability  of  Cylindrical  Shells  Under  Step  Loading" 

Points  emphasized: 

1.  Effect  of  axial  inertia  is  approximated. 

2.  Static  buckling  behavior  is  studied  for  the  discrete  model  used  in  the 
dynamic  analysis. 

3.  Experimental  results  are  used  for  values  of  the  initial  imperfections. 
Equations :  Nonlinear  Donnell's  equation  in  terms  of  displacement. 

Procedure: 

1.  Assume  appropriate  radial  displacement  function  w  and  initial  imperfection 
function  w. 

2.  Find  displacements  u  and  v  from  in-plane  equilibrium  equation  and  in 

turn  o  ,  o,  and  x 

x  y  xy 

3.  Evaluate  K.E.  &  P.E.  of  system. 

4.  Apply  Hamilton's  principle  resulting  in  a  nonlinear  ordinary  differential 
equation  (4  degrees  of  freedom). 

Static  Analysis  -  To  insure  existence  of  postbuckling  equilibrium  positions. 

In  Numerical  Study:  Parameters  to  be  chosen  are: 

1.  3  wave  numbers  i,  k,  £,  (i=2k) 

2.  Initial  imperfections. 

Dynamic  Stability: 

1.  Use  4  nonlinear  D.E.  and  initial  imperfections  under  a  special  step  load. 

2.  Use  a  Runge-Kutta  numerical  integration  scheme  for  the  starting  solutions. 

3.  Switch  to  a  6th  order  predictor-corrector  method. 

4.  Integrate  over  interval  to  obtain  desired  accuracy. 
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Results:  Dynamic  buckling  load  parameters  are: 

1.  Modes  of  radial  deformation  (according  to  wave  number  i,  k,  and  £). 

2.  Initial  imperfections. 

R  R 

3.  Geometric  parameters  and  material  constants  and  v. 

Axial  inertia  results  in  a  dynamic  buckling  load  much  lower  than  static  and  there¬ 
fore  must  be  included  in  the  model. 

"For  finite  length  shells,  reflections  of  the  axial  stress  wave  must  be  taken 
into  account",  (p.  193) 


Nondimensional  time  defined  as  T 
Critique: 


i  JT 

Rip 


1.  First  to  assess  the  influence  of  a  realistic  triggering  mechanism  on  the 
calculated  dynamic  buckling  loads. 

2.  First  to  include  in-plane  inertia. 

3.  Nonlinear  theory  used. 

4.  Solved  nonlinear  partial  differential  equations  and  examined  stability 
with  respect  to  wave  numbers  resulting  in  Mathieu  type  equations. 

5.  No  stiffeners  -  this  analysis  for  unstiffened  circular  cylindrical  shells 
only. 


B-37 


NSWC  TR  79-447 


Ref.  No.  32 


EFIMOV,  A.  B.,  et.  al. 


June  1974 


"Loss  of  Stability  of  a  Cylindrical  Shell  on  Longitudinal  Impact" 


This  brief  paper  describes  the  phenomena  of  stability  loss  of  a  cylindrical 
shell  on  longitudinal  impact  against  the  end  by  an  infinite  mass. 


Stability  loss  at  the  initial  instant  of  bulging  is  axisymmetric  and  is  analyzed 
using  linearized  Donnell-type  equations.  A  random  function  w  (x,  y)  is  used  to 
describe  initial  shell  irregularities.  The  equation  used  is: 


W  (x,  y,  t) 


wq(x,  y) 


where :  D 


Eh 


12(1  -  V  ) 


R  *  radius 
h  *  thickness 
v  «  poissons  ratio 


This  equation  is  solved  to  find  that  there  are  two  axisynsnetric  waves  traveling 
down  the  cylinder  of  which  one  becomes  unbounded  as  impact  velocity  increases 
beyond  a  critical  velocity  Vg. 


when  V  *  0.30  the  author  claims  that  experiments  have  shown  this  to  be  true  to 
within  12%. 

(No  Experimentation  presented  by  author  or  claimed  to  have  been  accomplished  by 
Author) . 
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1.  Uniform  impact  across  end. 

2.  Inelastic  collision. 

3.  Stiffeners  are  not  accounted  for. 

4.  Indicates  a  new  approach  to  the  stability  of  shells  on  longitudinal  impact. 
Note:  Mr.  Efimov  was  a  protege  of  Mr.  Grigolyulc,  Ref.  33. 
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Ref.  No .  33 


GRIGOLYUK,  E.  I. 


December  1975 


"Effect  of  Shape  of  External  Pressure  Impulse  and  Design  Parameters  of  Thin  Circular 
Cylindrical  Shells  On  Their  Stability" 


To  ignore  inertia  terms ,  the  rate  of  load  application  is  assumed  greater  than 
the  rate  of  deflection  growth. 

Bends  are  on  the  order  of  shell  thickness.  The  shell  is  viewed  as  an  elastic 
system  with  K  degrees  of  freedom. 

Illustrated  is  the  behavior  of  a  shell. 

Critique: 

1.  Avoids  inertia  terms. 

2.  No  stiffeners. 

3.  Purely  analytical  but  no  equations  or  programs  presented  in  the  paper. 
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Ref.  No .  34 


NSWC  TR  79-447 


KALNINS, 


( 


A. 


!  "Dynamic  Buckling  of  Axisymmetric  Shells" 

i 


December  1974 


The  object  of  this  report  is  to  extend  the  linear  dynamic  buckling  analysis 
for  arbitrary  axisymmetric  shells  when  subject  to  axisymmetric  periodic  loads 
of  long  durations. 

It  is  limited  to  cases  where  the  applied  load  is  of  sufficiently  long  duration 
when  compared  to  the  longest  resonant  period  of  the  shell. 

Starting  with  the  equation  of  motion  of  a  superimposed  disturbance  on  an 
initial  stressed  shell,  a  series  solution  for  the  displacement  is  substituted 
and  the  Galerkin  method  used  to  examine  flexural  modes  which  lead  to  approxima¬ 
tions  of  the  governing  equations  of  the  form  of  the  Hill  equation  and  result  in 
the  Mathieu  equation. 

Damping  is  incorporated  in  the  analysis  of  the  response  modes  of  the  shell 
by  incorporating  it  into  the  Mathieu  equation. 

The  case  of  the  Nuclear  reactor  cooling  pipe  shroud  is  examined. 


Critique: 

1.  Axisymmetric  shells  only. 

2.  Limited  dynamics. 

3.  Damping  included. 

4.  Linear  theory. 
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Ref.  No.  35 

GOROSHKO,  0.  A.  and  EMEL'YANENKO,  V.  V.  July  1975 


"Dynamic  Stability  of  Layered  Anisotropic  Shells” 


Linear  theory  is  used  to  study  the  effects  of  the  structure  of  the  material 
on  the  dynamic  stability  of  orthotropic  cylindrical  shells  subject  to  periodic 
loads  for  various  boundary  conditions. 

This  is  a  mathematical  treatise  filled  with  assumptions  to  make  higher  order 
terms  disappear.  At  different  times  both  Love's  and  Donnell's  theories  are  used 
resulting  in: 

Results :  Mathieu's  equation  and  three  regional  equations  for  a  given  set  of 
parameters,  a  set  of  curves  was  generated  showing  minimum  critical  loads. 

Critique: 

1.  Interesting  but  for  wrong  type  of  loading. 

2.  Otherwise,  the  result  seems  to  be  the  type  of  thing  we  are  looking  for. 

3.  Linear  theory. 
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Ref.  No.  36 

MAYMON  and  LIBAI  November  1977 


"Dynamics  and  Failure  of  Cylindrical  Shells  Subjected  to  Axial  Impact" 


Purpose :  1)  To  present  a  mathematical  model  of  the  behavior  of  cylindrical  shells 

under  axial  impact  including  statistical  characteristics  of  initial  imperfections. 
2)  To  present  a  digital  computer  program  parametric  study  of  different  shells. 

Method:  Linearized  Donnell-type  equations  derived  by  Baruch  and  Singer; 

-  closely  spaced  stiffeners  to  calculate  displacements  are  smeared  over 
the  shell  surface. 

-  at  least  2  stiffeners  must  be  in  each  wavelength. 

-  imperfections  approximated  by  Fourier  series  expansions. 

-  stress  calculations  -  discreteness  effects  taken  into  account, 

membrane  stresses 
bending  stresses 

-  normal  distribution  of  initial  imperfections  assumed. 

Failure  Criteria:  It  was  assumed  that  a  stress  oriented  design  limitation  would 
be  reached  before  any  significant  nonlinear  effects  developed. 

Results :  A  computer  program  capable  of  using  excessive  displacements  as  failure 
criteria.  The  complexity  of  nonlinear  analysis  of  stiffened  shells  is  so  great 
that  its  routine  use  is  not  practical  at  the  present  time. 

Axial  imperfections  are  most  critical  for  dynamic  buckling  than  the  circum¬ 
ferential  ones. 

More  effort  is  needed  in  establishing  real  statistical  imperfection  distribu¬ 
tions  and  experimental  verification  of  analytical  modes  of  dynamic  collapse. 

Critique: 


1.  Linear  theory  (see  Failure  Criteria) 

2.  Smeared  stiffeners. 

3.  This  is  the  most  comprehensive  paper  to  date. 

4.  Computer  program  very  large. 
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Ref.  No.  37 

LIBAI,  A,  and  DURBAN,  D.  December  1977 


"Buckling  of  Cylindrical  Shells  Subjected  to  Nonuniform  Axial  Loads" 


This  paper  deals  with  linear  buckling  of  a  cylindrical  shell  with  circumferen¬ 
tially  varying  axial  edge  loads  or  thermal  loads. 

Closed  form  expressions  for  the  eigenvalues  are  obtained. 

The  results  are  checked  against  a  numerical  study.  The  primary  concern  is 
with  heating  a  cylinder  along  a  strip  to  initiate  buckling. 

Linearized  Donnell  equations  used  in  conjunction  with  membrane  theory. 

Critique:  Not  remotely  near  what  we  are  looking  for,  but  is  an  example  of  the 
closed  form  solution  to  the  linearized  Donnell  equations. 
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Ref.  No.  38 


ZIMCIK,  D.  G. 


March  1976 


"Stability  of  Circular  Cylindrical  Shells  Under  Transient  Axial  Impulsive  Loading" 


This  paper  reports  the  results  of  experimental  testing  of  thin-walled  circular 
cylindrical  shells  subjected  to  dynamic  transient  axial,  square  wave  loading  of 
varying  time  duration. 

Types  of  shells  tested: 

-  geometrically  near  perfect 

-  axisymmetric  imperfect 

-  asymmetric  imperfect  (p~0.1) 

Special  manufacturing  techniques  were  used.  Material  was  birefringent  liquid 
epoxy. 

Results: 

1.  Indicated  an  increase  in  the  dynamic  buckling  stress  above  the  static 
value  for  all  shells  subjected  to  transient  square-pulse  loading. 

2.  Curves  were  drawn  showing  a  definite  finite-time  buckling  impulse 
minimum  in  the  time  domain  investigated. 

Shell  response  observed  experimentally  was  used  to  verify  predictions  obtained 
from  an  analytic  model: 

-  Karman-Donnell  compatibility  equilibrium  equations. 

-  Galerkin  procedures  used  to  get  set  of  nonlinear  partial  D.E. 

-  Modified  linear  stress  wave  analysis  used. 

Experimentation  and  analysis  show  good  agreement. 

Conclusions : 

1.  Buckling  stiffness  increases  for  short  duration  of  loading. 

2.  The  simplified  analytical  model  adequately  described  the  shells  buckling 
response. 

"The  two  quantities,  dynamic  buckling  load  and  pulse  duration 
can  be  combined  as  a  finite-time  buckling  impulse  to  provide 
a  conservative  dynamic  buckling  design  criterion",  (p.  25) 

This  was  first  study  providing  experimental  data  on  controlled  asymmetric 
shape  imperfection  effects  on  static  or  dynamic  buckling. 


B-45 


jimhLaaiy 


NSWC  TR  79-447 

Critique: 

1.  Model  testing  may  be  a  cost  effective  way  to  design  shells  to  meet 
dynamic  buckling  requirements. 

2.  Boundary  condition  effects  ignored  in  the  analysis. 

3.  No  stiffeners. 

4.  Inertia  included  (axial,  rotary,  and  radial). 

3.  Both  linear  and  nonlinear  theories  used. 

6.  Initial  imperfection  study  is  most  comprehensive  and  relative  to  date. 
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TERMS 

Rate  of  propagation  of  longitudinal  waves  in  the  shell, 

(E/ipU-v2)))5* 

Coefficients  used  in  imperfection  approximations 

Eh^/ 12( 1-v2 ) ;  flexural  rigidity 

Modulus  of  elasticity 

Airy  Stress  Function 

Shear  modulus 

Shell  wall  thickness 

Wave  number 

Buckling  stress  coefficient 
Shell  length 

Membrane  stress  resultants 
Critical  axial  load 
Axial  load 

Axial  load  coefficients 
Loads 

Shear  force  resultant 
Radius  of  shell  curvature 

Shell  thickness,  time,  or  nondimens ional  time  depending  on  the  author 
Initial  radial  shape  imperfection  function 

Displacements  in  the  axial,  circumferential,  and  radial  directions 
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TERMS  (Cont.) 

x,y,z  Axial,  circumferential,  and  radial  coordinates  respectively 

v  End  shortening  speed 

Z  Batdorf  Parameter,  nondimens ional  having  several  different  forms 

L2/Rt  or  I^/Rtd-V2)** 

v  Poisson's  ratio 

Axial  half  wavelength 

A  Circumferential  half  wavelength 

y 

V  Differential  operator 

3^/dx^  +  2(3<!*/3x23y2)  +  3^/3y^ 

a  Stress 

0  Critical  stress 

c 

6  Material  density,  generally 
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Appendix  C 

RESULTS  OF  COSMIC  COMPUTER  PROGRAM  SEARCH 
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Geometrically  Nonlinear  Static  and  Dynamic  Analysis  of  Arbitrarily  Loaded 

Shells  of  Revolution  (SADAALSOR) 

(Naval  Postgraduate  School) 


The  design  of  many  shell  structures  is  influenced  by  the  geometrically  nonlinear 
response  of  the  shell  when  subjected  to  static  and/or  dynamic  loads.  As  a  conse¬ 
quence,  a  number  of  investigations  have  been  devoted  to  the  study  of  the  buckling 
phenomenon  exhibited  by  shells.  Most  early  works  examine  the  behavior  of  the 
shallow  spherical  cap,  the  truncated  cone,  and  the  cylinder  under  axisymmetric 
loads.  Due  to  the  lack  of  information  on  the  axisymmetric  response  of  shells 
with  other  meridional  geometries  and  on  the  response  of  shells  subjected  to 
asymmetric  loads,  a  computer  program  for  the  geometrically  nonlinear  static  and 
dynamic  response  of  arbitrarily  loaded  shells  of  revolution  has  been  developed. 

The  program  can  be  used  to  analyze  any  shell  of  revolution  for  which  the 
following  conditions  hold: 

1.  The  geometric  and  material  properties  of  the  shell  are  axisymmetric,  but 
may  vary  along  the  shell  meridian. 

2.  The  applied  pressure  and  temperature  distributions  and  initial  conditions 
are  symmetric  about  a  datum  meridional  plane. 

3.  The  shell  material  is  isotropic,  but  the  modulus  of  elasticity  may  vary 
through  the  thickness.  Poisson's  ratio  is  constant. 

4.  The  boundaries  of  the  shell  may  be  closed,  free,  fixed,  or  elastically 
restrained. 

The  governing  partial  differential  equations  are  based  upon  Sanders'  nonlinear 
thin  shell  theory  for  the  condition  of  small  strains  and  moderately  small  rotations. 

At  each  load  or  time  step,  an  estimate  of  the  solution  is  obtained  by  extrapola¬ 
tion  from  the  solutions  at  the  previous  load  or  time  steps.  The  sets  of  algebraic 
equations  are  repeatedly  solved  using  Potter's  form  of  Gaussian  elimination,  and 
the  pseudo  loads  are  recomputed,  until  the  solution  converges. 

An  automatic  variable  load  incrementing  routine  is  included  in  the  program 
for  the  static  analysis.  Post-buckling  behavior  cannot  be  determined  in  the  static 
analysis  because  of  the  method  of  solution  employed. 
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Dynamic  Response  of  Shells  of  Revolution  During  Vertical  Impact  Into 
Water  -  Hydroelastic  Interaction  (SORVERT) 

(North  American  Aviation,  Inc.) 


This  is  one  of  several  programs  developed  out  of  a  study  to  determine  the 
hydro-elastic  response  of  representation  of  the  structure  of  the  Apollo  Command 
Module  immediately  following  impact  on  water. 

The  program  determines  the  hydroelastic  response  of  a  flexible  shell  of  revolu¬ 
tion  during  an  axially  symmetric  impact  into  an  incompressible  fluid,  and  is  based 
on  the  following  theory. 

During  the  impact,  the  deformation  of  the  body  (in  particular,  its  rate  of 
deformation)  interacts  with  the  fluid  flow  to  produce  oscillating  pressure  distri¬ 
butions  on  the  wetted  surface  of  the  body.  The  hydroelastic  response  of  the 
flexible  shell  is  obtained  by  the  numerical  solution  of  the  combined  hydrodynamic 
and  shell  equations.  The  results  obtained  are  compared  numerically  with  those 
derived  neglecting  the  interaction  and  applying  rigid-body  pressures  to  the  same 
elastic  shell. 
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HOLBOAT  -  General  Instability  Analysis  of  Inhomogeneous,  Anisotropic, 
Stiffened  Cylinders  Under  Combined  Loads 
(Martin  Marietta  Corp.) 


The  HOLBOAT  computer  program  provides  an  instability  analysis  of  inhomogeneous, 
anisotropic,  right  circular  cylinder  (or  segment)  under  combined  loading.  It 
is  an  efficient,  closed-form  solution  based  on  the  Kirchof f-Love  hypothesis, 
general  anisotropic  constitutive  equations,  and  Flugge's  differential  equations 
of  equilibrium.  Proper  use  of  this  program  is  based  on  the  user  being  familiar 
with  composites  technology  nomenclature  and  theory,  shell  instability  theory, 
and  FORTRAN  programming.  The  analysis  for  axial  loads,  pressure,  and  torsion 
is  based  upon  the  formulation  of  Cheng  and  Ho,  who  developed  a  solution  in  terms 
of  classical  small  deflection  theory  and  Flugge's  differential  equation  of 
equilibrium.  Their  analysis  has  been  extended  to  include  bending. 

Three  operating  modes  are  available  to  the  user: 

1.  Input  of  experimentally  determined  values  of  stiffness  from  isolated 
axial  compression,  torsion,  and  internal  pressure  test. 

2.  Direct  input  of  stiffness  matrices. 

3.  Input  of  cylinder  geometry  and  wall  construction. 

This  program  can  be  used  to  analyze  right  circular  cylinders  or  segments 
of  right  circular  cylinders.  Vertical  and  circumferential  stiffening  members 
may  be  included  in  the  modeling.  Stringers  or  rings  may  be  on  the  inside  or 
outside  of  the  cylinder.  Wall  construction  may  consist  of  a  skin  alone  or  a 
skin  reinforced  with  stiffening  members.  The  skin  may  be  a  laminated  cylinder 
constructed  of  layers  of  different  materials  having  different  elastic  properties 
or  orientations.  Each  individual  layer  may  be  isotropic  or  orthotropic;  symmetric 
or  balanced  arrangement  of  layers  is  not  required. 
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Appendix  D 

INITIAL  IMPERFECTIONS  SUMMARY 


Previous  methods  of  dealing  with  initial  imperfections: 

1.  A  four-parameter  modal  approximation 

u>  =  (H/R)  (dj  cosITm  x  cos  ny  +  d ^  cos  2mTTx)  K  =  xR/L  (29,30) 

0)  ■  h  (dj  cos  cos  +  d2  cos  -p-)  (12) 

x  y  x 

2.  Experimental  Results  -  actual  measurements  (24,31) 

or  constructed  to  give  specific  values  (38) 

3.  A  random  function  u)Q(x,y)  (32) 

4.  Approximated  using  Fourier  Series  Expansions  (36) 


U) 

o 


^  ^  ]  a  sin  -s—  cos  (-^  -  0  ) 

JLmJ  m,n  L  R  n 

m  n 


where  0^  =  phase  angle. 

5.  Asymmetric  initial  imperfections  assumed  (27) 

6.  Initial  imperfections  symmetric  about  a  plane  through  the  axia  are  accounted 
for  (16) 


7.  Approximated  by  "white  noise"  (13) 

8.  Imperfection  is  assumed  in  the  shape  of  the  buckling  mode  (3) 

9.  Assumed  form  of  initial  displacement  u>:  (2) 

e>  =  sin  (p-)  sin  (^-) 

Ax  Ay 

Notes:  The  numbers  given  in  parenthesis  after  each  method  refer  to  references 
given  in  the  body  of  the  report. 

Nomenclature  used  on  this  page  is  defined  in  Appendix  B. 
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